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quantum magnets, and the disordered contact process. 
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1. Preamble 

Phase transitions belong to the most fascinating phenomena in nature. The large 
scale structure of the universe is the result of a series of phase transitions during 
the early stages of its development. At low temperatures, quantum phase transitions 
between different ground states lead to unconventional behavior and exotic new phases 
of matter. Phase transitions also occur between nonequilibrium states, e.g., during 
growth processes or chemical reactions. Even our everyday life is unimaginable without 
the never ending transformations of water between ice, liquid and vapor. 

Under normal conditions, the phase transitions of water are so-called first-order 
transitions. At these transitions, the two phases coexist at the transition point, and 
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a finite amount of heat (the latent heat) is released when going from one phase to the 
other. Transitions that do not involve phase coexistence and latent heat arc called 
continuous transitions or critical points. They are particularly interesting because the 
typical length and time scales of fluctuations diverge when approaching the transition 
point. Understanding the resulting singidarities of physical observablcs, the so-called 
critical behavior has been a prime endeavor in physics, and the concepts established 
during this process, viz., scaling and the renormalization group, now belong to the 
central paradigms of modern physics. 

Many realistic systems contain a certain amount of quenched (i.e., time- 
independent or frozen in) disorder. This disorder can take the form of vacancies 
or impurity atoms in a crystal lattice, or it can consist of extended defects such as 
dislocations or grain boundaries. The question of how quenched disorder influences 
phase transitions and critical points is both conceptually very interesting and of 
fundamental importance for practical applications. It is therefore not surprising that 
this question has attracted considerable attention over the course of the last thirty 
years or so. Early work often focused on the average behavior of the disorder on larger 
and larger length scales. 

However, it has become increasingly clear that an important role is played by rare 
but large disorder fluctuations and the rare spatial regions where they occur. These 
regions can lead to a singularity (the so-called Griffiths singularity) in the free energy 
in an entire parameter region close to the phase transition that is now known as the 
Griffiths region or Griffiths phase. In generic classical systems, the thermodynamic 
Griffiths singularity is very weak and probably unobservable in experiment. In recent 
years, rare region effects have reattracted considerable attention, in particular in 
connection with nonequilibrium phase transitions and with zero-temperature quantum 
phase transitions. At these transitions, rare regions can have much stronger effects, 
ranging from from strong power-law singularities in the free energy to a complete 
destruction of the phase transition. 

The goal of this review is to provide a unifying framework for rare region effects 
on classical, quantum, and nonequilibrium phase transitions with quenched disorder. 
The core of the discussion will be limited to (order-disorder) phase transitions between 
conventional phases, i.e., we assume that the bulk phases are not changed qualitatively 
by the presence of the disorder. The microscopic effect of the impurities and defects 
then consists in locally favoring one of the two phases over the other. This type of 
quenched disorder is sometimes referred to as weak disorder, or random- Tc disorder 
(because it changes the local critical temperature at a thermal transition), or, from 
analogy with quantum field theory, random-mass disorder. 

The structure of this review is as follows. In section 2, we collect the basic concepts 
of phase transitions and critical phenomena necessary for the purpose of this review. 
In section 3 we give an introduction to the influence of quenched disorder on phase 
transitions with particular emphasis on rare region effects. We then put forward a 
general classification of weakly disordered phase transitions according to the effective 
dimensionality of the rare regions. In the bulk of the review, consisting of sections 4, 
5, and 6, we illustrate these general ideas by discussing in detail rare regions effects 
at several classical, quantum, and nonequilibrium phase transitions, respectively. We 
conclude in section 7 with a summary and a discussion of experiments and open 
questions. 
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2. Phase transitions and critical behavior 

In this section, we briefly collect the basic theoretical concepts of phase transitions and 
critical phenomena to the extent necessary for the purpose of this review. More details 
can be found in textbooks on this topic such as those by Ma [1] or Goldenfeld [2]. 

2.1. Landau theory 

Most modern theories of phase transitions are built on the foundation of Landau 
theory [3 6]. Landau introduced the general concept of an order parameter, a 
thermodynamic quantity that vanishes in one phase (the disordered phase|) and is 
non-zero in the other phase (the ordered phase). Very often the choice of an order 
parameter for a particular transition is obvious, e.g.. for the ferromagnetic transition 
where the total magnetization is an order parameter. Sometimes, however, finding an 
appropriate order parameter is a complicated problem in itself, e.g., for the disorder- 
driven localization-delocalization transition of non-interacting electrons. 

Landau theory can be understood as the unification of earlier mean-field theories 
such as the van-der-Waals theory of the liquid gas transition [7] or Weiss' molecular 
field theory of ferromagnetism [8]. It is based on the crucial assumption that the free 
energy is an analytic function of the order parameter m and can thus be expanded in 
a power series, 

F = FL{m) = Fo + rm.'^ + vmP + um.'^ + 0{m^) . (1) 

Here r,v,u are parameters that depend on all degrees of freedom other then the order 
parameter m. The physical value of m is the one that minimizes Fl- For sufficiently 
large r, the minimum free energy is always located at to = (the system is in the 
disordered phase), while for sufficiently small r, the minimum as at some m ^ 
(ordered phase). If u 7^ 0, the transition from to = to to ^ occurs discontinuously, 
i.e., it is of first order. If w = (as is often the case by symmetry), the theory describes 
a continuous transition or critical point at r = 0. Thus, r is measuring the distance 
from the critical point, r cx (T — Tc) for a thermal transition. 

Within Landau theory, the qualitative behavior of all critical points is identical. 
Specifically, the order parameter vanishes as to = (— r/2u)^/^ when the critical point 
is approached from the ordered phase. This is an example of the so-called super- 
universality of Landau theory: The critical exponent /3 which describes the singularity 
of the order parameter at the critical point via to cx |r|'^, is predicted to have the 
mean-field value 1/2 for all critical points. Universality of the critical exponents 
is actually observed in experiments but it is weaker than the super-universality 
predicted by Landau theory, and the exponent values are in general different from what 
Landau theory predicts. Moreover, the exponent values turn out to be dimensionality 
dependent. For instance, all three-dimensional Ising ferromagnets have a common 
P w 0.32 while two-dimensional Ising magnets have /3 = 1/8. Three-dimensional 
ferromagnets with 0(3) Heisenberg symmetry also have a common value oi (3 k, 0.35 
but it is different from the one in Ising magnets. 

The reason for the failure of Landau theory to correctly describe the critical 
behavior is that it does not adequately include the fluctuations of the order parameter 
about its average value. The effects of these fluctuations in general decrease with 

X Here "disordered phase" refers to a phase without long-range order. This is not to be confused 
with quenched disorder in the form of impurities and defects. 
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increasing dimensionality and with increasing number of order parameter components. 
This c'xplains why the critical behavior of Ising magnets (one order parameter 
component) deviates more strongly from Landau theory than that of Heisenberg 
magnets (three components), and why the three-dimensional exponents are closer 
to the mean-field values than those of two-dimensional systems. 

This suggests that Landau theory might actually be correct for systems with 
sufficiently high dimensionality d. In fact, the fluctuations lead to two different critical 
dimensionalities, rf^ and dj for a given phase transition. If the dimensionality d is 
larger than the upper critical dimension fluctuations are unimportant for the 
leading critical behavior, and Landau theory gives the correct answers. If d is between 
the upper and the lower critical dimensions, rf^ > > d", a phase transition still 
exists but the critical behavior is different from mean-field theory. For dimensionalities 
below the lower critical dimension, fluctuations become so strong that they completely 
destroy the ordered phase. For the ferromagnetic transition at nonzero temperature, 
d+ = 4, and = 2 or 1 for Heisenberg and Ising symmetries, respectively. 

2.2. Scaling and the renormalization group 

As indicated in the last subsection, below the upper critical dimension order 
parameter fluctuations play a crucial role in determining the critical behavior. To 
include them, one has to generalize the Landau free energy (1) by writing the partition 
function Z as a functional integral 

Z = e-^l^ = J D[<^]e-^[^1 , (2) 

with the action or Landau-Ginzburg- Wilson (LGW) functional S given by 

^[<^] = ^j d''x [c(V<^(x))2 + FL{ct>{K)) - hcl>{x)] . (3) 

(j) is a fluctuating field whose average value is equal to the order parameter, m = {4>). 
(Here, (. . .) denotes the average with respect to the statistical weight e~^^'^^). We have 
also included an external field h conjugate to the order parameter. 

In the disordered phase, away from the critical point, the correlation function of 
the order parameter fiuctuations, G(x — y) = (0(x)(/)(y)) is generically short-ranged. 
When the critical point is approached, the correlations become long-ranged. Their 
typical length scale, the correlation length ^, diverges when the distance r from the 
critical point vanishes§, 

C^\rr ■ (4) 

Here, v is the correlation length critical exponent. Close to the critical point, the 
correlation length ^ is the only relevant length scale in the system. Therefore, the 
physical properties must be unchanged, if we rescale all Ic^ngths in the system by a 
common factor 6, and at the same time adjust the external parameters in such a way 
that the correlation length retains its old value. This gives rise to the homogeneity 
relation for the free energy density / = — (T/V) InZ, 

f{r,h) = b-'^f{rb^''',hby'^). (5) 

§ One needs to distinguish between the bare value of r which appears in the LGW functional (3) 
and the jjliysical or renormalized value which measures the distance from the true critical point. We 
will explicitly make this distinction when necessary but otherwise suppress it. 
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Table 1. Definitions of tiic commonly used critical exponents, m = (0) is the 
order parameter, and h is tfie conjugate field, r denotes the distance from the 
critical point and d is the space dimensionality. (The exponent yii defined in (5) 
is related to 5 by 2//1 = dS/{l + S).) 



The scale factor b is an arbitrary positive number, and ph is another critical exponent. 
Analogous homogeneity relations for other thermodynamic quantities can be obtained 
by differentiating /. These homogeneity laws were first obtained phenomenologically 
by Widom [9] and are sometimes summarily called the scaling hypothesis. Within the 
framework of the modern renormalization group theory of phase transitions [10] the 
scaling laws can be derived from first principles. 

In addition to the critical exponents v and defined above, a number of other 
exponents is in common use. They describe the dependence of the order parameter 
and its correlations on the distance from the critical point and on the field conjugate 
to the order parameter. The definitions of the most commonly used critical exponents 
are summarized in Table 1. Note that not all the exponents defined in Table 1 are 
independent from each other. The four thermodynamic exponents a, (3, 7, 5 can all 
be obtained from the free energy (5) which contains only two independent exponents. 
They are therefore connected by the so-called scaling relations 

2 - a = 2/3 + 7 , (6) 

2-a = /3((5 + l) . (7) 

Analogously, the exponents of the correlation length and correlation function are 
connected by the relations 

2-a = dv, (8) 

7 ={2-Ti)v . (9) 

Note that scaling relations explicitly involving the dimensionality, the so called 
hyperscaling relations, only hold below the upper critical dimension c?+. Their 
breakdown above is caused by dangerously irrelevant variables, see, e.g., Ref. [2]. 

In addition to the diverging length scale ^, a critical point is characterized by a 
diverging time scale, the correlation time ^4. It leads to the phenomenon of critical 
slowing down, i.e., very slow relaxation towards equilibrium near a critical point. At 
generic critical points, the divergence of the correlation time follows a power law 
oc where z is the dynamical critical exponent. At some transitions, in particular 
in the presence of quenched disorder, the divergence can be exponential, In^t oc 
The latter case is referred to as activated dynamical scaling in contrast to the generic 
power-law dynamical scaling. 
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As mentioned above, the critical exponents display the remarkable phenomenon 
of universality, i.e., they arc the same for entire classes of phase transitions which 
may occur in very different physical systems. These classes, the so-called universality 
classes, are determined only by the symmetries of the Hamiltonian and the spatial 
dimensionality of the system. The physical mechanism behind universality is the 
divergence of the correlation length. Close to the critical point the system effectively 
averages over large volumes rendering the microscopic details of the Hamiltonian 
unimportant. 

2.3. Classical versus quantum phase transitions 

In classical statistical mechanics, statics and dynamics decouple. This can be seen 
by considering a classical Hamiltonian H{pi,qi) = i?kin(Pi) + Hpot{qi) consisting of a 
kinetic part iJkin that only depends on the momenta Pi and a potential part iJpot that 
only depends on the coordinates qi. The classical partition function of such a system, 

^ = / n dp^e-"'''"/^ J n dq^e-^p-/^ = Zkin^pot , (10) 

factorizes in kinetic and potential parts which are completely independent from 
each other. The kinetic contribution to the free energy density will usually not 
display any singularities, since it derives from a product of simple Gaussian integrals. 
Therefore one can study the thermodynamic critical behavior in classical systems using 
time-independent theories like the Landau-Ginzburg- Wilson theory discussed in the 
last subsection. As a result of this decoupling, the dynamical critical exponent is 
completely independent from the thermodynamic ones. 

In quantum statistical mechanics, the situation is different. The Hamiltonian 
H and its parts iJkin and iJpot are now operators, and i/kin and i/pot in general do 
not commute. Consequently, the partition function Z — Tre^^^^ does not factorize, 
and one must solve for the dynamics together with the thermodynamics. Quantum 
mechanical analogs of the Landau-Ginzburg- Wilson theory (3) therefore need to be 
formulated in terms of space and time dependent fields; they can be derived, e.g., via 
a functional integral representation of the partition function. The simplest example 
of such a quantum Landau-Ginzburg- Wilson functional takes the form 

S[<^]= j^J^ dT J d^x [a(a,0(x,r))2 + c(V0(x,r))2 + Fi(<^(x,r))-/7,<^(x,T)] (11) 

with r being the imaginary time. This functional also illustrates another remarkable 
feature of quantum statistical mechanics. The imaginary time variable r effectively 
acts as an additional coordinate. For nonzero temperatures, 1/T < oo, the extra 
dimension extends only over a finite interval. If one is sufficiently close to the critical 
point so that the condition > 1/T is fulfilled, the extra dimension will not affect 
the leading critical behavior. We thus conclude that the asymptotic critical behavior 
at any transition with a nonzero critical temperature is purely classical. However, 
a quantum phase transition at T = is described by a theory that effectively is 
in a different (higher) dimension. Therefore, it will be; in a different universality 
class. At finite temperatures and in the vicinity of a quantum phase transition, the 
interplay between classical and quantum fluctuations leads to nontrivial crossover 
phenomena [11-14]. 
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Using the equivalence between imaginary time and an additional dimension, the 

homogeneity law for the free energy can be easily generalized to the quantum case 
(see, e.g., [12]). For conventional power law dynamical scaling it takes the form 

f{r,h,T) = b-^'^+''\f{rb^/'',hby'',Tb'') . (12) 

This relation implies that the scaling behavior of a zero-temperature quantum phase 
transition in d dimensions can be mapped onto that of some classical transition in 
d + z spatial dimensions. If space and time enter the theory symmetrically the 
dynamical exponent is z = 1, but in general, it can take any positive value. Note 
that this quantum to classical mapping is valid for the thermodynamics only. Other 
properties like the real time dynamics at finite temperatures require more careful 
considerations [12]. Moreover, some quantum phase transitions, such as the metal- 
insulator transition, do not have a classical counterpart and thus the quantum to 
classical mapping does not directly apply. 

Let us also emphasize that the entire concept of a Landau-Ginzburg- Wilson theory 
in terms of the order parameter field relies on the order parameter fluctuations being 
the only soft or gapless mode in the system. If there are other soft modes, e.g., 
due to conservation laws or broken continuous symmetries, they lead to long-range 
power-law correlations of various quantities even away from the critical point. This 
phenomenon is called generic scale invariance [15-17]. If one insists on deriving a 
Landau-Ginzburg- Wilson theory for a phase transition in the presence of generic scale 
invariance, the additional soft modes lead to coeSicients that are singular functions 
of space and time, severely limiting the usefulness of the theory [18-20]. In such a 
situation one should instead work with a coupled theory that keeps all soft modes at 
the same footing. This question has been explored in detail in Ref. [14]. In practise, 
complications due to generic scale invariance are more common for quantum phase 
transitions than for classical ones because (i) there are more soft modes at T = 
than at r > 0, and (ii) phenomena that classically only affect the dynamics influence 
the static critical behavior of a quantum phase transition as well. In addition to 
this mechanism, there are other reasons for a quantum phase transition being more 
complicated than one might naively expect, including topological effects. In many 
of these cases, a description of the transition in terms of a Landau-Ginzburg- Wilson 
theory is not possible. 

2.4- Finite-size scaling 

We now turn to the question of how a finite system size influences the behavior at a 
critical point, a question that will be very important later on for our discussion of rare 
region effects. In general, a sharp phase transition can only exist in the thermodynamic 
limit, i.e., in an infinite system. A finite system size results in a rounding and shifting of 
the critical singularities. If the extension of the system is finite only in some directions 
and infinite in others, a sharp phase transition can still exist, but the critical point will 
be shifted compared to the bulk value. Moreover, the universality class will change. 

The behavior of finite size systems in the vicinity of a critical point is described by 
finite-size scaling theory [21-23]. Finite-size scaling starts from the observation that 
the inverse system size acts as an additional parameter that takes the system away 
from the critical point. Because the correlation length of the infinite system £^ao is the 
only relevant length scale close to the critical point, finite-size effects in a system of 
linear size L must be controlled by the ratio L/^oo only. We can therefore generalize 
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the classical homogeneity relation (5) for the free energy density by including the 
system size 

/(r, h, L) = h-^fir 6l/^ /i i ft-^). (13) 
The quantum version (12) of the homogeneity law can be generalized in the same way: 

f{r,h,T,L) = b-^''+''>f{rb^/'',hby\Tb'-,Lb-^) . (14) 

We note in passing that the temperature scaling in the quantum homogeneity law 
can be viewed as finite-size scaling in imaginary time direction with the inverse 
temperature playing the role of the system size. 

The shift of the critical temperature or critical coupling as a function of L for 
geometries that still allow a sharp phase transition at finite L follows directly from 
(13). Setting b = L, h = 0, we obtain 

f{r,L)=L''F{rL'/n (15) 

where F{x) is a dimensionless scaling function. The quantum case can be treated 
analogously. The finite-L phase transition corresponds to a singularity in the scaling 
function at some nonzero argument Xc- The transition thus occurs at rcL^^" = Xc, 
and the transition temperature Tc{L) of the finite-size system is shifted from the bulk 
value T° by 

T,{L)-T° ^re = XeL-^^'' . (16) 

Note that the simple form of finite-size scaling summarized above is only valid 
below the upper critical dimension (i+ of the phase transition. Finite-size sc;aling can 
be generalized to dimensions above d+ , but this requires taking dangerously irrelevant 
variables into account. One important consequence is that the shift of the critical 
temperature, Tc{L) — T'c ~ L~'^ , is controlled by an exponent (j) which in general is 
different from 

3. Quenched disorder effects 

In this section, we give a brief introduction into the effects of quenched disorder on 
phase transitions and critical points. We focus on what is arguably the simplest and 
most benign type of disorder, viz., impurities and defects that lead to spatial variations 
of the coupling strength (i.e., spatial variations in the tendency towards the ordered 
phase) but not to frustration or random external fields. This also implies that the two 
bulk phases (that are separated by the transition) are not changed qualitatively by 
the presence of the quenched disorder. This type of quenched disorder is sometimes 
referred to as weak disorder, random- Tc disorder, or, in analogy to quantum field 
theory, random-mass disorder. 

In ferromagnetic materials, random- Tc disorder can be achieved, e.g., by diluting 
the lattice, i.e., by replacing magnetic atoms with nonmagnetic ones. Within a 
Landau- Ginzburg- Wilson theory such as (3) or (11), random-Tc disorder can be 
modelled by making the bare distance from the critical point a random function 
of spatial position, r — s- r -|- 5r{x). In this representation, the character of the 
impurities and defects is encoded in the statistical properties of the random-Tc term 
5r(x). Point-like impurities correspond to short-range correlations of (5r(x) in all 
directions while linear and planar defects lead to perfect correlations of (5r(x) in 
one and two dimensions, respectively. Let us also briefly comment on the disorder 
probability distribution P{6r). As long as the physics is dominated by long- wavelength 
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properties and the average behavior of the disorder, the details of the probabiUty 
distribution should not play an important role. In theoretical investigations a realistic 
distribution is often replaced by a Gaussian because the latter can be easily handled 
mathematically. We emphasize, however, that in some cases the dominating physics 
occurs in the tail of the probability distribution. In these cases, the form of the tail is 
important, and more careful considerations are required. 

Adding weak, random- Tc, quenched disorder to a system undergoing a continuous 
phase transition naturally leads to the following questions: 

• Will the phase transition remain sharp in the presence of quenched disorder? 

• If so, will the critical behavior change quantitatively (different universality class 
with new exponents) or even qualitatively (exotic non-power-law scaling)? 

• Will only the transition itself be influenced or also the behavior in its vicinity? 
3.1. Average disorder and Harris criterion 

The question of how quenched disorder influences phase transitions has a long 
history. Initially it was suspected that disorder destroys any critical point because 
in the presence of defects, the system divides itself up into spatial regions which 
independently undergo the phase transition at different temperatures (see Ref. [24] 
and references therein). However, subsequently it became clear that generically a 
phase transition remains sharp in the presence of defects, at least for classical systems 
with short-range disorder correlations. 

Harris [25] derived a criterion for the perturbative stability of a clean critical 
point against weak disorder. Harris considered the effective local critical temperature 
in blocks of linear size ^ which is given by an average of r + 5r{x) over the volume 
y = He observed that a sharp phase transition can only occur if the variation Ar of 
these local critical temperatures from block to block is smaller than the global distance 
from the critical point r. For short-range correlated disorder, the central limit theorem 
yields Ar oc ^~'^/^ oc r'^"/'^. Thus, a clean critical point is perturbatively stable, if the 
clean critic;al cixponents fulflU the inequality r'^'^^'^ < r for r — > 0. This implies the 
exponent inequality 

di^>2 (17) 

which is called the Harris criterion. Note that the Harris criterion is a necessary 
condition for the stability of a clean fixed point, not a suflicient one. It only deals 
with the average behavior of the disorder at large length scales. However, effects due 
to qualitatively new physics at finite length scales (and finite disorder strength) are 
not covered by the Harris criterion. 

The Harris criterion can be used as the basis for a classification of critical points 
with quenched disorder, based on the behavior of the average disorder strength with 
increasing length scale, i.e., under coarse graining. Three classes of critical points can 
be distinguished [26]. 

(i) The first class contains systems whose clean critical points fulfill the Harris 
criterion du > 2. At these phase transitions, the disorder strength decreases under 
coarse graining, and the system becomes asymptotically homogeneous at large length 
scales. Consequently, the quenched disorder is asymptotically unimportant, and 
the critical behavior of the dirty system is identical to that of the clean system. 
Technically, this means the disorder is renormalization group irrelevant at the clean 
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critical fixed point. In this class of systems, the macroscopic observables are self- 
averaging at the critical point, i.e., the relative width of their probability distributions 
vanishes in the thermodynamic limit [27, 28] . A prototypical example in this class 
is the three-dimensional classical Heisenberg model whose clean correlation length 
exponent is » 0.698 (see, e.g., [29]), fulfilling the Harris criterion. 

If a clean critical point violates the Harris criterion, it is destabilized by weak 
quenched disorder, and the behavior must change. Nonetheless, a sharp critical point 
can still exist in the presence of the disorder, but it must fall into one of the two 
following classes. 

(ii) In the second class, the system remains inhomogeneous at all length scales 
with the relative strength of the inhomogcncitics approaching a finite value for large 
length scales. The resulting critical point still displays conventional power-law scaling 
but with new critical exponents which differ from those of the clean system (and fulfill 
the inequality dv > 2). These transitions are controlled by rcnormalization group 
fixed points with a nonzero value of the disorder strength. Macroscopic observables 
are not self-averaging, but in the thermodynamic limit, the relative width of their 
probability distributions approaches a size-independent constant [27,28]. An example 
in this class is the classical three-dimensional Ising model. Its clean correlation length 
exponent, u « 0.627 (see, e.g. [30]) does not fulfill the Harris criterion. Introduction 
of quenched disorder, e.g., via dilution, thus leads to a new critical point with an 
exponent of ~ 0.684 [31]. 

(iii) At critical points in the third class, the relative magnitude of the 
inhomogencities increases without limit under coarse graining. The corresponding 
rcnormalization group fixed points are characterized by infinite disorder strength. 
At these infinite-randomness critical points, the power-law scaling is replaced by 
activated (exponential) scaling. The probability distributions of macroscopic variables 
become very broad (even on a logarithmic scale) with the width diverging with system 
size. Consequently, their averages are often dominated by rare events, e.g., spatial 
regions with atypical disorder configurations. This type of behavior was first found 
in the McCoy- Wu model, a two-dimensional Ising model with bond disorder perfectly 
correlated in one dimension [32,33]. However, it was fully understood only when 
Fisher [34, 35] solved the one-dimensional random transverse field Ising model by a 
version of the Ma-Dasgupta-Hu real space rcnormalization group [36]. Since then, 
several infinite-randomness critical points have been identified, mainly at quantum 
phase transitions since the disorder, being perfectly correlated in (imaginary) time, 
has a stronger effect on quantum phase transitions than on thermal ones. Examples 
include one-dimensional random quantum spin chains as well as one-dimensional and 
two-dimensional random quantum Ising models [26,37 40]. Very recently, a similar 
infinite-randomness fixed point has been found in the disordered one-dimensional 
contact process [41-43]. 

3.2. Rare regions and Griffiths effects 

In the last subsection we have discussed scaling scenarios for phase transitions with 
quenched disorder based on the global, i.e., average, behavior of the disorder strength 
under coarse graining. In this subsection, we turn to the main topic of this review, 
viz., the effects of rare strong spatial disorder fluctuations. 

Let us start by considering the example of a diluted classical ferromagnet as 
sketched in figure 1. The dilution reduces the tendency towards magnetic long-range 
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Figure 1. Sketch of a diluted magnet. The shaded region is devoid of impurities 

and therefore acts as a piece of the clean bulk system 



order and thus reduces Tc from its clean bulk value T^. However, in an infinite 
system one can find arbitrarily large spatial regions that are devoid of impurities. For 
temperatures between T° and Tc, these regions will show local magnetic order even 
though the bulk system is globally in the paramagnetic phase. These spatial regions 
are known as rare regions and the order parameter fluctuations induced by them belong 
to a class of excitations known as local moments; sometimes they are also referred to 
as instantons. The dynamics of the rare regions is very slow because flipping them 
requires a coherent change of the order parameter over a large volume. Griffiths [44] 
was the first to show that these rare regions can lead to a singularity in the free energy, 
the Griffiths singularity, in the entire temperature region Tc < T < Tc which is now 
known as the Griffiths region or the Griffiths phase [45] . Analogous singularities also 
exist on the ordered side of the phase transition. 

The probability w for finding a rare region is exponentially small in its volume 
Vrr and in the impurity concentration p; up to pre-exponential factors it is given 
by w ~ exp(— pVrr). Rare regions are thus nonperturbative degrees of freedom that 
are not accounted for in conventional approaches to phase transitions and critical 
points that are based on perturbation theory and the perturbative renormalization 
group. They can be viewed as examples of a broader class of rare event phenomena 
that range from the well known Lifshitz tails in the density of states of disordered 
semiconductors [46,47] to the slowing down of computer algorithms in combinatorial 
optimization [48]. 

In classical systems with uncorrelated or short-range correlated disorder, 
thermodynamic Griffiths effects are very weak because the singularity in the free 
energy is only an essential one [49-51]. To the best of our knowledge, classical 
thermodynamic Griffiths singularities have therefore not been verified in experiment 
(see also [52]). In contrast, the long-time dynamics inside the Griffiths phase is 
dominated by the rare regions, as is signified by a non-exponential time-dependence 
of the spin autocorrelation function [45,53-56]. 

Long-range disorder correlations can increase the rare region effects qualitatively. 
In particular, if the disorder is perfectly correlated in some spatial directions, or - in 
the case of a zero-temperature quantum phase transition - in imaginary time, the rare 
regions are extended objects which are infinite in the correlated space or imaginary 
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time directions. This makes their dynamics even slower and so increases their effects. 
This enhancement of rare region effects was first observed in the abovcmcntioncd 
McCoy- Wu model [32,33], a two-dimensional classical Ising model with linear defects. 
More recently, it has been studied in great detail in various disordered quantum 
systems: In onc-dimcnsional and two-dimensional random transverse field Ising 
models, the singularity of the free energy in the Griffiths phase takes a power-law 
form with nonuniversal continuously varying exponents. Several thermodynamic 
observablcs including the average susceptibility actually diverge in a finite region of 
the disordered phase [26,34,35,39,40] rather than only at the critical point. Similar 
phenomena have also been found in quantum Ising spin glasses [57-59]. In addition, 
power-law Griffiths effects, sometimes also called Griffiths-McCoy effects, were also 
found at the nonequilibrium phase transition in a disordered one-dimensional contact 
process [41-43]. 

Recently, it has been shown that rare region effects can become strong enough 
to completely destroy the sharp phase transition. This happens, if the dynamics of 
the rare regions stops completely, i.e., if a static order parameter develops on an 
isolated rare region independent of the bulk system. The global phase transition is 
then smeared or rounded, as was found in itinerant quantum magnets [60] where 
the overdamping of the magnetic excitations by the electronic degrees of freedom 
suppresses the dynamics of sufficiently large rare regions [61,62]. Later, analogous 
smeared transitions been found in a variety of systems, ranging from classical Ising 
magnets with planar defects [63, 64] to nonequilibrium spreading transitions in the 
contact process [65,66]. 

3.3. A classification of rare region effects 

In the last subsection, we have seen that the influence of rare regions on continuous 
phase transitions can range from exponentially weak free energy contributions to a 
complete destruction of the sharp phase transition by smearing. It would clearly 
be very desirable to determine what conditions lead to which type of rare region 
effects. Recently, a general classification of rare region effects at order-disorder phase 
transitions in systems with weak, random-mass disorder has been suggested [67]. This 
classification is based on the following idea. The probability w for finding a large 
spatial region of volume Vrr ^ L^p devoid of impurities decreases exponentially 
with its volume, «;(Lrr) ~ exp(— pVrr). The importance of large rare regions thus 
depends on how rapidly the contribution of a single region to observable quantities 
increases with its size. It turns out that this question is controlled by the effective 
dimensionality ^rr of the rare regions, or more precisely, by the relation between cJrr 
and the lower critical dimension d~ of the phase transition under consideration, jj 

To understand this, let us consider a single isolated region of linear size Lrr 
that is locally in the ordered phase, i.e., the bulk distance from the critical point is 
negative, r — r(oo) < 0. We now ask: How does the renormalized distance from the 
critical point, r(LRR), (or equivalently, the energy gap) of the finite size rare region 
depend on its linear size Lrr? Three cases can be distinguished (a summary is given 
in table 2). 

II When counting the effective dimensionahties for quantum phase transitions we include the 
imaginary time direction as one of the dimensions. 
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Class 


RR dimension 


Griffiths effects 


Dirty critical point 


Critical scaling 


A 


C^RR < d~ 


weak cixponcntial 


conventional 


power law 


B 


dRR = d~ 


strong power-law 


infinite randomness 


activated 


C 


dRR > d~ 


RR become static 


smeared transition 


no scaling 



Table 2. Classification of rare region (RR) effects at critical points in the presence 
of weak quenched disorder according to the efltective dimensionality ^rr of the 
rare regions. For details see text. 



3.3.1. Class A: If the effective dimensionality dRR of the rare regions is below the 
lower critical dimensionality d~ of the problem, an isolated rare region cannot undergo 
the phase transition by itself. The renormalized r(LRR) is therefore positive, but it 
decreases with Lrr following a power law. The leading contributions of a rare region 
to thermodynamic quantities are controlled by r(LRR), e.g., the order parameter 
susceptibility x proportional to r^^{LRR). The contributions of a rare region to 
observables can therefore at most grow as a power of its linear size Lrr. (In the case 
of a classical magnet, % ~ S^g ~ ^rr)- T^^^s power-law increase cannot overcome the 
exponential drop in the rare region density w(Lrr). 

As a result, for transitions in this class, thermodynamic rare region effects are 
exponentially weak and characterized by an essential singularity in the free energy 
[44,49 -51]. Right at the dirty critical point, the rare region effects arc subleading, and 
the critical behavior is of conventional power-law type. Examples for critical points 
in class A can be found in generic classical equilibrium systems with point defects 
(where the rare regions are finite in all directions and thus c^rr = 0) but also at some 
quantum phase transitions such as the transition in the diluted bilayer Heisenberg 
quantum antiferromagnet [68-70] (here dRR = 1 but d~ = 2 so that the condition for 
class A is fulfilled). 

3.3.2. Class B: In this class, the rare regions are exactly at the lower critical 
dimension, rfRR = d~ . They still cannot undergo the phase transition independently, 
but the renormalized distance from the critical point (or equivalently, the energy gap) 
decreases exponentially with the rare region volume Vrr. Therefore, the rare region 
contribution to observables such as the susceptibility increases exponentially with size 
and can potentially overcome the exponential decrease of the rare region density. 

Generically, this results in a power-law low-energy density of states leading 
to strong power-law Griffiths singularities with a nonuniversal continuously varying 
exponents. The scaling behavior at the dirty critical point itself is dominated by the 
rare regions resulting in exotic activated (exponential) scaling instead of conventional 
power-law scaling. Class B is realized, e.g., in a classical Ising model with linear 
defects, i.e., disorder perfectly correlated in one direction (McCoy- Wu model [32,33]), 
in percolation with linear defects [71], and in random quantum Ising models (where 
the disorder correlations are in imaginary time direction) [26,34,35,39,40]. In all these 
systems, ^rr = d~ = 1. Several thermodynamic observables including the average 
susceptibility actually diverge in a finite region of the disordered phase rather than just 
at the critical point. Class B rare region effects also occur at the nonequilibrium phase 
transition of the contact process with point defects. Including the time dimension, 
this transition also has dRR = d~=l. 
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3.3.3. Class C: Finally, in class C, the rare regions can undergo the phase transition 
independently from the bulk system, i.e., they arc effectively above the lower critical 
dimension of the problem, > d~ . In this case, the dynamics of the locally ordered 
rare regions completely freezes, and they develop a truly static order parameter. As 
a result, the global phase transition is destroyed by smearing [60] because different 
spatial parts of the system order at different values of the control parameter. In 
the tail of the smeared transition, the order parameter is extremely inhomogeneous, 
with statically ordered islands or droplets coexisting with the disordered bulk of the 
system. This behavior was identified in itinerant quantum Ising magnets [60-62] where 
the damping due to the electronic degrees of freedom leads to an effective long-range 
(1/r^) interaction in time direction. At this transition d^iB. = 1 and d~ ~ \ but 
it nonetheless belongs to class C because the rare region can undergo a Kosterlitz- 
Thouless type phase transition to the ordered phase. Smeared transitions can also 
be fomid in classical magnets with planar defects [63,64] (c?rr = 2 > 1), in 

nonequilibrium spreading transitions with extended (linear or planar defects) [65, 66] 
where c^rr > 1 but d~ = 1 (including the time dimension) or in stratified percolation 



This classification is expected to apply to all continuous order-disorder transitions 
between conventional phases that can be described by Landau-Ginzburg- Wilson type 
theories with random-mass disorder and short-range interactions. It also relies on the 
assumption that the interaction between the rare regions can be neglected if their 
density is sufficiently low (i.e., sufficiently far away from the true disordered critical 
point, if any). While this seems a reasonable assumption in the case of short-range 
spatial interactions, long-range interactions will likely modify the rare region effects. 
In the case of itinerant quantum Heisenberg magnets, it has been argued that the 
interaction between the rare regions due to the long-range RKKY interaction strongly 
enhances the rare region effects, leading to a smeared transition even though a single 
rare region would not independently undergo the phase transition [73]. 

The classification of rare region effects provides a general framework for a variety 
of existing results, and it helps organizing the search for new phenomena. In the 
following sections of this review we will discuss in detail several examples of rare regions 
effects at classical, quantum, and nonequilibrium phase transitions with quenched 
disorder and relate them to the above classification. 

4. Classical phase transitions 

This section is devoted to rare region effects at classical equilibrium phase transitions 
in systems with random- Tc type disorder. We first consider the case of uncorrelated 
disorder (point-like defects) which leads to exponentially weak "classical" Griffiths 
effects. We then show that disorder correlation enhance the rare region effects by 
discussing systems with linear and planar defects. 

4..I. Randomly diluted Ising model 

For dcfiniteness, let us consider a randomly site-diluted classical Ising magnet on a 
d-dimensional hypercubic lattice oi N — L'^ sites. The Hamiltonian is given by 



[72] (rfRR = 2,d- = l). 




(18) 



{id) 
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Figure 2. Schematic phase diagram of a randomly diluted Ising model as function 
of temperature T and impurity concentration p. pc is the geometric percolation 
threshold of the lattice, and is the clean critical temperature 



Here, Si = ±1 is the Ising spin at lattice site i, h is the external magnetic field, J > 
is the exchange interaction, and indicates that the sum is over pairs of nearest 
neighbor sites only. The site dilution is introduced by the quenched random variables 
Ci with ei = 1 representing an occupied site and = corresponding to a vacancy. 
In our case of point defects, the are statistically independent with a probability 
density 

F(e) = (l-p)<5(e-l)+p5(e) , (19) 

i.e., p is the probability for finding a vacancy at a particular site. A schematic phase 
diagram of this system as function of p and temperature T is shown in figure 2. At 
zero temperature, magnetic long-range order survives for all concentrations p < Pc, 
where Pc is the geometric percolation threshold of the lattice. If p > Pc, no long-range 
order is possible as the system is decomposed into independent finite-size clusters. For 
p < Pc, magnetic order survives up to a nonzero temperature Tc(p); but at pc, it is 
destroyed immediately by thermal fluctuations^, Tc{pc) = [74-76]. 

In a seminal paper, Griffiths [44] showed that in the thermodynamic limit, the 
magnetization (per site), m — j^M — CjiSi, is a singular function of the magnetic 

field ft, at /i = for any temperature T below the transition temperature of the 
clean system (below the dashed line in figure 2) independent of p. Conventional 
paramagnetic behavior only appears above T^. This singularity and the corresponding 
singularities in other thermodynamic quantities are now known as the Griffiths 
singularities and the temperature region over which they occur is called the Griffiths 
region or the Griffiths phase [45]. Griffiths' proof of the singular dependence of the 
magnetization on the field relied on an application of the theorem of Yang and Lee on 
the complex zeros of the partition function [77,78], and it did not immediately give 
the functional form of the predicted singularity. However, after some controversy it 
was soon established that the Griffiths singularity is a weak essential singularity with 
the magnetization remaining infinitely differentiable dX h — Q [49,50,52]. 

^ For Ising symmetry, this appears counterintuitive as the fractal dimension of the critical percolation 
cluster is larger than the lower critical dimension = 1. However, for fractals, topological factors 
like ramification are important for magnetic ordering [74] . 
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4-2. Optimal fluctuation theory 

The leading singularities in the thermodynamic quantities in the Griffiths region can be 
easily estimated using optimal fluctuation theory arguments similar to those employed 
in the study of band tails of doped semiconductors [79,80]. 

In the site diluted Ising model (18), the probability w(Lrr) for finding a region 
of volume Lf^^^ devoid of vacancies is (to exponential accuracy) 

w(Lrr) ~ (1 - p)^^R = eM-pLRK) (20) 

with p = — ln(l — p). Let us now consider a temperature T in the paramagnetic 
Griffiths region, T^ip) <T <T°, ov equivalently ro = {T - T°)/T° < 0. The bare 
distance r(LRR) from the local critical point for a rare region of size Lrr is shifted by 
finite size cfFects because the rare region acts as a finite-size piece of the clean system. 
If ro is not too close to zero (in the Gaussian fluctuation region), the shift can be 
determined by finite size scaling (16) with the mean-field value v = 1/2, leading to 

r{Lnn)=ro+A/Li^ (21) 

where A is a constant. 

If r{Ljiji) < 0, the rare region is locally in the ordered phase. Thus, at a given 
To, only rare regions larger than L„iin(?'o) = ko/^|~^^^ will show local order. Since 
they are of finite size, the locally ordered rare regions do not undergo a true transition 
but coherently fluctuate due to thermal fluctuations. To estimate their contribution 
^RR to the free energy, we note that the contribution of a single region varies at most 
like a power of its size. To leading exponential accuracy, the free energy of the rare 
regions can therefore be estimated by summing their probabilities, 

POO 

Fkr ~ dLnR w{Lnn) ~ exp(-p|roM|-'^/2) . (22) 

Close to the bulk critical point ro — 0, i.e., in the clean bulk critical region, the 
mean fleld relation for the finite-size shift of r is not valid anymore and has to be 
replaced by finite size scaling with the shift exponent (j) = 1/v where u is the clean 
bulk correlation length exponent. The behavior of the free energy then crosses over 
toFRR~exp(-p|roM|-<*-). 

The response to a small external magnetic field can be derived in a similar way: 
In the Griffiths region, an isolated, locally ordered rare region essentially acts as a 
single big Ising spin of magnitude 5efi ~ ^rr- Thus its magnetization in an external 
field h is given by 

«(L„H) = S.,tanh(t5../T)^{"« "^^ZlHZ^ P3) 

where 6 is a constant that depends on tq. The singularity in the magnetization- field 

relation arises from the large rare regions whose magnetization is saturated by the 
field h. To exponential accuracy, the singular rare region contribution to the total 
magnetization is therefore given by 

/•oo 

Mrr= / dLRR w(Lrr) ~ exp(-p67/i) (24) 

in agreement with references [50, 52]. More recently, many of the properties of 
Griffiths phases have been proven rigorously in the mathematical physics literature, 
see e.g., [81-83] and references therein. 
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At this point, it is important to emphasize that the Griffiths singularities do 
not display the same degree of universality as the critical singularities right at a 
continuous phase transition. This is caused by the fact the rare regions are phenomena 
at finite length scales in the tail of the disorder probability distribution. Therefore, 
the functional form of the Griffiths singularities can depend on the specifics of the 
disorder distribution. We will now illustrate this non-universality by comparing the 
results (22) and (24) for the randomly diluted Ising model with those of a Landau- 
Ginzburg- Wilson theory with Gaussian disorder. 



4-3. Random-Tc Landau- Ginzhurg- Wilson theory 

In this subsection, we consider rare region effects in a classical Landau-Ginzburg- 
Wilson (LGW) theory with random- Tc disorder. From a critical phenomena point of 
view, the phase transition in this model belongs to the same universality class as the 
transition in the randomly site diluted Ising model considered above. Nonetheless, 
the functional form of the Griffiths singularities will turn out to be different from (22) 
and (24). 

Let us start from a d-dimensional LGW theory of the form [see equations (1) and 

(3)] 

S[d)] = [ d-^x [(V(/)(x))2 + (ro + (5r(x))(A2(x) + m,/>4(x) - /i(/>(x)] .(25) 



The disorder is modelled by the random contribution (5r(x) to the distance to the 
critical point. The distribution of fr(x) is usually chosen to be a Gaussian 



P[(5r(x)] (X exp 



j d'^x <5V(x) 



(26) 



In this model, a rare region is a large spatial region in which the typical value of 
ro + (5r(x) is negative while the bulk distance ro from the critical point is positive. 
One qualitative difference to the randomly diluted Ising model case is immediately 
obvious: Because the Gaussian disorder is unbounded, locally ordered rare regions 
can exist for any value of the bulk distance ro from the critical point. This means, 
there is no conventional paramagnetic phase, and the Griffiths region covers the entire 
high-temperature part of the phase diagram. 

The contribution of these rare regions to the free energy can be estimated 
using optimal fluctuation arguments [84] similar to those in the subsection 4.2. 
The probability for finding a large region of linear size Lrr and a typical value 
^Ttyp = — ro — s of the disorder variable is exponentially small, 



w(Lrr, s) ~ exp 



(27) 



Treating the finite-size corrections in mean-field theory (valid outside the asymptotic 
bulk critical region) as before, the bare distance from the critical point for such a rare 
region is given by 

r(LRR, s) = ro + Srtyp + A/Ll^ = -s + A/LI^ . (28) 

If r(LRR, s) < 0, the rare region is locally in the ordered phase. Thus, the minimum 
size permitting local order is -Lmin = (A/s)^/^. To leading exponential accuracy, 
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the contribution of the locally ordered rare regions to the free energy can again be 
estimated by summing their probabilities, 

PCO POO POO 

Frr~ / ds dLnRw{LnR,s)^ dsexp[-Cs-'^/^{s + rof] . (29) 

where C is a constant. The second equality arises because the integral over L is 
dominated by its lower bound. The integral over s can be estimated using the saddle 
point method. The leading contribution comes from ssp = ''oc'/(4 — d). This implies 
that with increasing ro, i.e., with increasing distance from the true transition, the 
dominating rare region contribution to the free energy comes from smaller and deeper 
disorder fluctuations. This is qualitatively different from the case of the randomly 
diluted Ising model where with increasing ro larger and larger islands dominate. 
Inserting the saddle point value into the equation for Frr gives 

Frr ~ expl-C'r^^-"^)/'] . (30) 

where C is another constant. This result has also been derived more rigorously by 
considering replica instanton solutions of the replicated and disorder averaged Landau- 
Ginzburg- Wilson theory [84,85]. 

Comparing the Griffiths singularities of the randomly diluted Ising model (22) 
and the LGW theory with Gaussian disorder (30), we find them both exponentially 
weak. However, the functional form of the singularity is not universal, it does depend 
on the specifics of the microscopic disorder distribution. Dotsenko [84] also calculated 
the rare region contribution to the free energy in a nonzero magnetic field for the 
LGW theory with Gaussian disorder. He found (5Frr(/i) ~ exp [-Eh-'^^^] with E a 
constant. To exponential accuracy, the same same singularity should also govern the 
magnetization. As in the zero-field case, the functional form of the Griffiths singularity 
is different from that of the randomly site diluted Ising model, equation (24). 

4-4- Dynamic rare region effects 

So far, we have focusscd on the effects of the rare regions on thermodynamic properties, 
and we have found them to be exponentially weak in generic classical systems with 
short-range correlated disorder. In contrast, rare region effects on the dynamics are 
much more dramatic, as will be discussed in this subsection. 

For definiteness, let us again consider the dilute ferromagnet (18). The dynamics 
is assumed to be purely relaxational, corresponding to model A in the classification of 
Hohenberg and Halperin [86]. Microscopically, this type of dynamics can be realized, 
e.g., via the Glauber [87] or Metropolis [88] algorithms. Because the rare region 
effects are local in space, their dynamic properties can be characterized in terms of 
the average spin autocorrelation function C{t) defined by 

c{t) = ^j2iSi{t)sm (31) 

i 

where Si{t) is the value of the spin at time t and (•) denotes the thermodynamic 
average. In the conventional paramagnetic phase T > (see figure 2), the average 
autocorrelation function decays exponentially fast because none of the rare regions is 
locally ordered. In contrast, in the Griffiths region Tc{p) < T < T°, we expect the 
rare regions to display slow dynamics. Following arguments arguments first made by 
Dhar [53] and later refined by Randeria et al. [45] and Bray [56,89], the rare region 



Rare region effects at classical, quantum, and nonequilibrium phase transitions 20 



part of the autocorrelation function can be estimated as the sum of independent 
contributions from each of the rare regions 

POO 

CRR(t)~ / d^RR w(iRR) exp[-V6(iRR,r-o)] . (32) 

Here w(iRR) is the probability for finding a rare region of size Lrr as given in (20) 
and Ct(-^RRj''o) is the relaxation time of such a rare region at reduced temperature 
ro- Depending on ro, equation (32) predicts different long time behavior of the 
autocorrelation function. 

(i) In the conventional paramagnetic phase, T > (tq > 0), the relaxation time 
remains finite even for the largest islands. Close to it behaves as (I/rr, ro) ~ ~ 
rQ^'^ for sufficiently large islands (Lrr > £, )■ Here ^ is the clean bulk correlation length 
and z and ly are the clean bulk critical exponents. The asymptotic time dependence 
of C{t) in this regime is a simple exponential. 

InC(f) - -t/r ■ (33) 

(h) At the clean critical temperature T = (ro = 0), i.e., at the boundary of 
the Griffiths region, the relaxation time of the large rare regions is given by finite- 
size scaling [22]. It diverges as ^((Lrr) ^ -Z^rr- After inserting this into (32), the 
integral over Lrr can be estimated using the saddle-point method. The asymptotic 
time dependence of C{t) is a stretched exponential 

\nC{t) ~ fd/{d+z) _ (34^ 

(iii) Inside the Griffiths region, Tc{p) < T <T^, the relaxation time is limited by 
the time taken to coherently reverse the entire rare region. Such a reversal requires 
creating a domain wall whose free energy is aL^^. Therefore, the relaxation time 
of an island of linear size Lrr is given by S,t{LRR) ~ Toexp[aL'^]. The factor 1/T 
in the Arrhenius behavior has been absorbed into the surface tension a, and tq is a 
microscopic time scale. For large t, the integral (32) can be calculated using the saddle 
point method leading to an even slower time dependence of the form 

InC(t) ~ -(lni)''/(''-i) . (35) 

More recently, this asymptotic decay of the autocorrelation function was established 
rigorously for disordered discrete lattice spin systems with bounded disorder 
distributions [90,91]. It should be noted however, that in Monte-Carlo simulations 
[92 95] the time dependence of C{t) seems to follow a stretched exponential form 
InC(t) ~ —t^ rather than the asymptotic form (35). This suggests that the asymptotic 
regime is only reached after a very long transient time interval. We also point out 
that (35) is universally valid only for bounded disorder distributions; unbounded 
distributions can lead to asymptotic forms different from (35) as was shown in Ref. [96] . 

So far our discussion of dynamic rare region effects has focused on Ising systems. 
For continuous spin symmetry, the behavior of the autocorrelation function at and 
above the clean critical temperature t'^ remains to be given by relations (33) and 
(34). However, in the Griffiths region, Tc{p) < T < T°, the behavior changes. The 
relaxation time of a rare region is much shorter than for Ising spins since there is no 
free energy barrier hindering relaxation. Rather, the relaxation occurs by diffusion 
of the order parameter, driven by thermal noise. Following Bray [56], the relaxation 
time can be estimated as follows: The change of cluster magnetization during some 
small time interval St is 5M. ~ L'^^ since the noise at different sites adds incoherently. 
Since several of such steps also act incoherently, the change of magnetization after 
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time t is SM.(t) ~ t^/^L^'^. Complete relaxation has occurred for d'M.{t) ~ M ~ 
giving £_t ~ ^RR- Putting this into (32) and evaluating the integral by the saddle point 
method gives 

InC(t) - , (36) 

a much faster decay than the corresponding equation (35) for the Ising case. The 
same result has also been derived in the large- A'' limit (where A'' refers to the number 
of spin components) [55]. In contrast to the Ising case, Montc-Carlo simulations of a 
diluted Heisenberg magnet [93] are in good agreement with (36) suggesting that the 
asymptotic regime is reached much faster in the Heisenberg case than in the Ising 
case. 



4-5. Disorder correlations and enhanced rare region effects 

In addition to localized defects such as impurity atoms or vacancies, many realistic 
systems also contain extended defects such as dislocations or grain boundaries. They 
can be modelled by quenched disorder that is perfectly correlated in one (linear defects) 
or two (planar defects) dimensions. 

It is interesting to ask whether such disorder correlations increase or decrease the 
influence of the defects. Naively, one might expect that disorder which is perfectly 
correlated in rfcor dimensions but uncorrelated in the remaining rfran = d — dcov 
dimensions is weaker than uncorrelated disorder because the system is "random only 
in some of the directions" . However, this naive assumption turns out to be incorrect. 
In fact, the opposite is true: Long-range disorder correlations enhance the influence of 
the defects because "it is harder to average out an extended fluctuation than a local 
one". 

The first systematic study of a phase transition in the presence of extended 
defects was carried out more than 30 years ago by McCoy and Wu. In a scries 
of papers [32,33,97-100], these authors studied a disordered two-dimensional Ising 
model, now known as McCoy- Wu model, in which all horizontal bonds are identical 
while all vertical bonds in each row are the same. However, the vertical bonds vary 
from row to row as independent random variables. Thus, the disorder is perfectly 
correlated in one dimension (the horizontal direction) but uncorrelated in the other 
(the vertical direction). This implies dcor = c?ran = 1- 

The McCoy- Wu model in zero magnetic field is partially exactly solvable with 
transfer-matrix type methods (see also [101]). McCoy and Wu [33] analyzed the 
free energy and found that there is only an essential singularity at the critical 
temperature. No exact results could be obtained for the spontaneous magnetization, 
but McCoy [98,99] calculated the properties of the spins on the surface of a half- 
plane cut perpendicular to the random direction. Because of Griffiths' theorem [102], 
magnetization and susceptibility of these boundary spins provide lower bounds for the 
corresponding bulk values. McCoy found the surprising result that the susceptibility is 
infinite for an entire range of temperatures above the critical temperature Tc- He also 
found that there is an even larger range of temperatures in which the susceptibility 
exists but the magnetization is still a nonanalytic function of the field. These unusual 
results initially gave rise to some confusion as to whether the phase transition in the 
McCoy- Wu model is sharp or smeared. 

By now, it has become clear that the transition is actually sharp, i.e., the 
spontaneous magnetization develops as a collective effect of the whole system at the 
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critical temperature which is also the (only) temperature at which the spatial 
correlation length becomes infinite. The unusual behavior of the susceptibility in the 
vicinity of the critical temperature is a rare region effect that can be understood as 
follows: For simplicity, let us assume that the random vertical bonds in the McCoy- Wu 
model have a binary probability distribution 

P{J) = {l-p)S{J-Jo)+p 6{J - cJo) (37) 

with c < 1. Thus, p is the probability of a vertical bond being weak. In an infinite 
sample, one can find arbitrarily large "stacks" of rows with only strong vertical bonds. 
The probability for finding such a rare region of width Lrr is given by 

w{Lnn) ~ (1 - P)^^^ ~ exp(-pLRR) (38) 

with p = — ln(l — p). In the Griffiths region below the transition temperature 
of the clean system with only strong vertical bonds, these rare regions are locally in 
the ordered phase. Since each rare region is infinite in the horizontal direction and of 
finite width Lrr in the vertical direction, it is equivalent to a one-dimensional Ising 
model with an effective interaction proportional to Lrr. As such, it cannot undergo a 
phase transition by itself but its susceptibility increases exponentially with its width, 

x{Lrr) ~ ^ exp(aLRR) (39) 

where the constant a vanishes at the clean critical temperature and increases with 
decreasing temperature (see, e.g., [2]). 

The rare region contribution to the total susceptibility can be obtained by 
summing over all rare regions, 

Xnn-ldL^n.iL^n)x{Lnn). (40) 

This integral diverges when a > p before the true critical point of the dirty system is 
reached (the divergence does not require coherence between the rare regions). Thus, 
in the Griffiths region of the McCoy- Wu model, the exponential rarity of large rare 
regions is overcome by an exponential increase of a rare region's susceptibility with 
its size. The probability distribution Wy^ of local susceptibilities Xl can be obtained by 
combining (38) and (39), leading to 

^Mxi)-x.;'^''-' (41) 

where z' = a/p is an exponent that varies continuously throughout the Griffiths region. 
In contrast, the corresponding probability distribution for the the case of uncorrelated 
disorder decreases exponentially for large xi because the susceptibility of a rare region 
only increases like a power of Lrr, see equation (23). The perfect disorder correlations 
in the McCoy- Wu model therefore enhance the rare region effects from exponentially 
weak classical Griffiths singularities to much stronger power-law singularities that are 
sometimes called Griffiths-McCoy singularities. 

Let us emphasize that none of the peculiar properties of McCoy and Wu's exact 
solution are reflected in the standard perturbative renormalization group approach to 
the problem. Early renormalization group work [103] based on a single expansion 
in e = 4 — d did not produce a critical fixed point. This was interpreted as a 
smeared transition or a first order one [104,105]. An alternative approach, the double-e 
expansion in powers of e = 4— d and the number dcoi of correlated dimensions [106-108] 
gave rise to a critical fixed point with a conventional power-law singularity in the 
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free energy rather than the essential singularity found in McCoy and Wu's solution. 
However, the general arguments on the importance of rare regions made in section 
3.3 point to a qualitative difference between d^or = 1 in the McCoy- Wu model and 
rfcor 1 assumed in the double-e expansion. For rfcor = 1, the rare regions (whose 
effective dimension djiji = dcor) are at the lower critical dimension rf", leading to the 
exponential size dependence (39) of the susceptibility. In contrast, for dcor < 1, the 
rare regions are below d~. Thus, the susceptibility depends on the rare region size 
via a power law, and the Griffiths effects in the double-e expansion are exponentially 
small. Therefore it is not surprising that the double-e expansion approach cannot 
capture the physics of the McCoy- Wu model. 

The partition function of the McCoy- Wu model written in terms of the transfer 
matrix in horizontal direction is essentially equivalent to the partition function of a 
quantum mechanical random transverse-field Ising model [12,101]. Much insight can 
therefore be obtained by studying this quantum Ising model. In fact, the behavior of 
the McCoy- Wu model remained incompletely understood until Fisher [34, 35] applied 
a version of the Ma-Dasgupta-Hu [36] real-space renormalization group to the random 
transverse-field Ising model. We will discuss these developments in more detail in 
sections 5.1 to 5.3. 



4-6. Smeared phase transition in an Ising model with planar defects 

In the last subsections we have seen that perfect disorder correlations in one direction 
greatly enhance the rare region effects. In this subsection, we consider planar defects, 
i.e., disorder perfectly correlated in two dimensions. Recently, it has been shown 
[63,64] that they have even more dramatic consequences: They completely destroy 
the sharp Ising phase transition by smearing. 

Most of our arguments below will be rather general. If necessary for definiteness, 
we focus on a three-dimensional Ising model with nearest-neighbor interactions on a 
cubic lattice. In the clean system, all interactions are identical and have the value 
Jq. The defects are modelled via 'weak' bonds randomly distributed in one dimension 
(uncorrelated direction). The bonds in the remaining two dimensions (correlated 
directions) remain equal to Jq. The system effectively consists of blocks separated by 
parallel planes of weak bonds. Thus, d^an = 1 and dcor = 2. The Hamiltonian of the 
system is given by: 

H = — JiSij^kSi+ij,k — Jo{Si,j,kSij+i,k + Sij^kSij,k+i), (42) 

i=l,. . . , Lran *=!)■■■ i^ran 

j,k = l LcoT J)fc = l)---.i'cor 

where iran(-^cor) is the length in the uncorrelated (correlated) direction, i, j and k are 
integers counting the sites of the cubic lattice, and Ji is the random coupling constant 
in the uncorrelated direction. The J, are drawn from a binary distribution: 

P{J) = {l-p)S{J-Jo)+pS{J-cJo) (43) 

characterized by the concentration p and the relative strength c of the weak bonds 
(0 < c< 1). 

Let us now consider rare regions in this system which consist of "thick slabs" 
containing only strong bonds in the uncorrelated direction. As in the case of point 
or linear defects, for temperatures below the clean critical temperature T°, these rare 
regions are locally in the ordered phase even if bulk system is still in the disordered 
(paramagnetic) phase. However, the behavior of locally ordered planar rare regions 
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Figure 3. Schematic behavior of the magnetization in the tail of a smeared phase 
transition for bounded and unbounded disorder distributions. For details see text. 

differs qualitatively of that of localized or linear rare regions. Each planar rare region 
is infinite in the two correlated dimensions but finite in the uncorrelated direction. 
Thus, it is equivalent to a two dimensional fsing model that can undergo a real 
phase transition independently of the rest of the system. Thus, each rare region 
can independently develop true static order with a non-zero static value of the local 
magnetization. (Griffiths theorem [102] ensures that the interaction of the rare region 
with the paramagnetic bulk system can only increase its magnetization over that of a 
completely isolated "slab".) Once static order has developed, the magnetizations 
of different rare regions can be aligned by an infinitesimally small interaction or 
external field. The resulting phase transition will thus be markedly different from 
a conventional continuous phase transition. At a conventional transition, a non-zero 
order parameter develops as a collective effect of the entire system which is signified 
by a diverging correlation length of the order parameter fluctuations at the critical 
point. In contrast, in a system with planar defects, different parts of the system 
(in the uncorrelated direction) will order independently, at different temperatures. 
Therefore the global order will develop inhomogeneously and the correlation length 
in the uncorrelated direction will remain finite at all temperatures. This defines a 
smeared or rounded phase transition. 

Note that similar to the Griffiths effects discussed in the proceeding sections, 
the properties of the smeared transition depend on whether or not the disorder 
distribution is bounded (see figure 3). If the disorder distribution is bounded like 
the binary distribution (43), there is a true paramagnetic phase with zero total order 
parameter at high temperatures. At some temperature Trr, (identical to the clean 
critical temperature for the distribution (43)), a non-zero order parameter starts to 
develop on the rare regions, accompanied by an essential singularity in the free energy 
density (which stems from the probability for finding a rare region). In contrast, for 
distributions that are unbounded in the sense that they permit rare regions with a 
local Tc = oo (e.g., a Gaussian distribution), the total order parameter is non-zero for 
all temperatures, and the free energy density is analytic. 

4-7. Optimal fluctuation theory for smeared phase transition 



The leading thermodynamic behavior in the tail of the smeared transition can be 
determined using optimal fluctuation theory similar to section 4.2. Following [63], 
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we develop this theory for the general case of rfcor > 2 correlated dimensions 

and dran ^ d — dcoi uncorrclatcd dimensions. Wc first consider bounded disorder 
distributions like the binary distribution (43) and later briefly describe the differences 
for Gaussian disorder. 

The probability w of finding a thick slab of width Lrr devoid of any weak bonds 
is, up to pre-exponcntial factors, given by 

w(iR,R,) - cxpi-pL'^^") . (44) 

with p = — ln(l — p). Such a rare region develops static long-range order at some 
temperature TdL^i^) below the clean critical point T^. The value of Tc(Lrr) varies 
with the size of the region: The largest islands develop long-range order closest to the 
clean critical point. As discussed in section 2.4, finite size scaling for the clean system 
(because the rare regions are free of defects) yields 

TO - T,(Lrr) = |re(LRR)| = AL-^ (45) 

where (j) is the finite-size scaling shift exponent and A is the amplitude for the 
crossover from d dimensions to a slab geometry infinite in dcor dimensions but finite 
in dj-an = d — dcor dimensions. If the total dimensionality d is below the upper critical 
dimension d+ = 4, hypcrscaling is valid and the shift exponent </> = The reduced 
temperature r = T—T^ measures the distance from the clean critical point. Combining 
(44) and (45) we obtain the probability for finding a rare region which becomes critical 
at T*(; as 

w{rc) ~ exp(-B |rc|"*""/"^) (for r ^ 0-) (46) 

where tli(^ (;onstant B is given by B = pA'^'""^/'^. The total (or average) order 
parameter m is obtained by integrating over all rare regions having rc(iRR) > r. 
Since the functional dependence on r of the order parameter on a given island is of 
power-law type it only enters the pre-exponential factors. Therefore we obtain to 
exponential accuracy 

m(r) - exp(-B |r|-'^""/'^) (for r ^ 0-) . (47) 

Once a nonzero density of ordered and aligned rare regions exist, the Ising symmetry 
is spontaneously globally broken. The rare regions produce an effective background 
magnetic field everywhere in space which cuts off any possible further singularities. 

The spatial magnetization distribution in the tail of the smeared transition is 
very inhomogeneous. On the already ordered islands, the local order parameter m(x) 
is of the same order of magnitude as in the clean system. Away from these islands 
it decays exponentially with the distance from the nearest island. The typical local 
order parameter mtyp can be estimated via the typical distance of any point from the 
nearest ordered island. Prom (46) we obtain 

xtyp ~ exp(B Irl-'^-Z'^/dran) ■ (48) 

At this distance from an ordered island, the local order parameter has decayed to 

mtyp ~ e""^*^"/^" ~ exp 

where is the bulk correlation length (which is finite and changes slowly throughout 
the tail region of the smeared transition) and C is constant. A comparison with (47) 
gives the relation between mtyp and the thermodynamic order parameter m, 

|logmtyp| ~m~^/*»" . (50) 



-Cexp(B |r| 



(49) 
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Figure 4. Left; Logarithm of the magnetization m as a function of — T\~^ 
{v = 0.627) for impurity concentrations p = 0.2,0.25,0.3 and strengtli c = 0.1. 
Inset: Decay slope B as a function of — log(l — p). Riglit: Local magnetization mi 
as function of position i in the uncorrclated direction at temperature T = 4.425 
(one particular disorder realization). Lower panel: Coupling constant Ji in the 
uncorrelated direction as a function of position i. Inset; Log-linear plot of the 
region in the vicinity of the largest ordered island (from [64]). 

Thus, TOtyp decays exponentially with m indicating an extremely broad distribution 
of the local order parameter TOx = {Sy^). Its functional form has been determined in 
reference [63], it reads 

P[log(mx)] ~ I log(mx)|*"""^ (for rrix < 1) . (51) 

These predictions of the optimal fluctuation theory have been tested by large 
scale computer simulations. A model with infinite-range interactions in the correlated 
directions was studied in reference [63]. This mean- field type calculation permitted 
large system sizes of up to Lran = 10^ making it very suitable for identifying effects of 
rare events. The results were in excellent agreement with the theoretical predictions. 
Note that a similar mean-field model had been investigated earlier [109] but for much 
smaller sizes of up to Lran = 1000. The raw data in this paper are very similar to 
reference [63] , btit they were interpreted in terms of power-law critical behavior with 
an unusually large order parameter exponent fi « 3.6. This may bo partially due to 
an unfortunate choice of parameters (in particular, a higher impurity concentration of 
p = 0.5) which makes it hard to extract the functional form of the magnetization tail. 

Recently, the short-range Ising model with planar defects (42) has been studied 
by large scale Monto-Carlo simulations [64] using the Wolff cluster algorithm [110]. 
The results confirm the smeared transition scenario and the predictions of optimal 
fluctuation theory for a realistic model. A typical result for the magnetization in the 
tail of the smeared transition below » 4.511 is presented in the left panel of figure 
4. It shows the logarithm of the total magnetization vs. ]T° — T\~'' averaged over 
240 samples for system size Z/ran = 200, icor = 280 and three disorder concentrations 
p = {0.2.0.25,0.3}. Here, u = 0.627 is the three-dimensional clean critical Ising 
exponent. For all three concentrations the data follow (47) over more than an order 
of magnitude in m. The deviation from the straight line for small m is due to the 



Rare region effects at classical, quantum, and nonequilibrium phase transitions 27 



conventional finite size effects. In the inset we show that the decay constant B depends 
linearly on p = — log(l — p). 

The right panel of figure 4 illustrates the very inhomogeneous character of the 
local magnetization in the tail of the smeared transition. It displays the local 
magnetization nii of a particular disorder realization as a function of the position 
i in the uncorrelated direction for the size Lran = 200, Lcor = 200 at a temperature 
T = 4.425 in the tail of the smeared transition. The lower panel shows the local 
coupling constant Ji as a function of i. The figure shows that a sizable magnetization 
has developed on the longest island only (around position i = 160). One can also 
observe that order starts to emerge on the next longest island located close to i = 25. 
Far form these islands the system is still in its disordered phase. 

As discussed in the last subsection, the functional form of the tail of the smeared 
transition is nonuniversal. In particular, for unbounded disorder distributions, it 
differs from (47). In the case of Gaussian disorder, the magnetization for large r 
behaves as m ~ exp(— Sr^"'*™-/'^), i.e., he magnetization tail reaches all the way to 
T = DO [63]. 

4-8. Dynamics at a smeared phase transition 

The static magnetic order on the rare regions in an Ising model with planar defects 
also modifies the dynamics. The leading long-time behavior can be determined by 
generalizing the approach of section 4.4 to the case of extended defects, as was done 
in reference [111]. The behavior above and at the clean critical temperature T° turns 
out to be identical to the conventional dynamic Griffiths effects discussed in section 
4.4. This is not surprising because above there is no difference between the Griffiths 
and the smearing scenarios: All rare regions are locally still in the disordered phase. 
Specifically, the asymptotic time dependence of the spin autocorrelation function C{t) 
right at the clean critical temperature is given by 

InC(t) p^/(dra„+2) ^dra„/(d«„+z) (52) 

where z is the clean critical exponent of the d-dimensional system. For temperatures 
below T|?, the behavior at a smeared transition differs from the dynamic Griffiths 
effects of section 4.4. Here, at short times before a crossover time tx, the decay of 
C(t) is still given by by the stretched exponential (52). For t > tx, the system realizes 
that some of the rare regions have developed static order and contribute to a non-zero 
equilibrium value C(oo) of the autocorrelation function C{t). The approach of C{t) 
to this equilibrium value follows a nonuniversal power law 

C(t) - C(oo) - t-'f' . (53) 

The value of tp could not be found since it depends on the pre-exponential factors 
neglected in the optimal fluctuation theory. 

Reference [111] also discussed the rare region contribution to the nonequilibrium 
relaxation of the magnetization in the tail of the smeared transition. It turned 
out to be controlled by decay laws with the same functional form as those of the 
autocorrelation function. Those results were confirmed by computer simulations of a 
infinite-range model similar to that in reference [63]. Figure 5 shows a typical result 
of these simulations. The left panel shows the time evolution of the magnetization of 
a large system of Lj-an = 10''' sites, averaged over 25 disorder realizations, for several 
ttmipciratures from T = 2.8 to 3.25 in the vicinity of the clean critical point which 
is r" = 3 for the infinite-range model. The behavior at T° follows the stretched 
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Figure 5. Left: Time evolution of tlie magnetization for several temperatures 
from T = 2.8 to 3.25 (top to bottom) in the vicinity of the clean critical point 

= 3. The impurity concentration is p = 0.5 their strength is c = 0.2. The 
runs start from a fully polarized state at ta = and go to tmax = 2000. Right; 
Double-logarithmic plot of the approach of the magnetization to the equilibrium 
value for T = 2.90, 2.95 and 2.98 (top to bottom). The straight lines are fits to 
power laws yielding exponents of 1.78, 2.20, and 3.31, respectively (from [111]). 



exponential him{t) ~ — t^/^ as expected from the analog of (52) because the mean- 
field value of the clean dynamical exponent is z = 2, and dran = 1- For T > T^^, 
the magnetization decay is faster, and a detailed analysis [111] shows it to be simple 
exponential. For T < Tj^, the magnetization approaches a nonzero equilibrium value 
TOeg. The approach to this equilibrium value is analyzed in the right panel of figure 
5. The difference m{t) — m^q follows a nonuniversal power law m{t) — m^q ^ t~'^ . 

5. Quantum phase transitions 

In section 4, we have seen that rare region effects in classical systems with generic 
short-range-correlated disorder are exponentially weak and that stronger effects 
require long-range disorder correlations. In this section, we discuss quantum phase 
transitions which occur at zero temperature and thus require the inclusion of the 
imaginary time direction, as discussed in section 2. Quenched impurities are time 
independent, the disorder is therefore always correlated in at least one of the relevant 
dimensions, and we generically expect rare region effects at a quantum phase transition 
to be stronger than at a classical transition. 

The main focus of this section will again be on order-disorder transitions between 
conventional phases that can be described by Landau-Ginzburg- Wilson theories. 
However, the concept of quantum rare region effects or quantum Griffiths singularities 
is much broader, and in section 5.7, we will briefly discuss other examples. 

5.1. Random transverse field Ising model 

One of the prototypical models (if not the prototypical model) displaying a quantum 
phase transition is the transverse field Ising model. Its Hamiltonian 

ff = -^J.,5f5|-^/ir5f (54) 
ihj) » 
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consists of two noncommuting parts, the first being the exchange interaction between 
the z-componcnts of the spin-^ operators Si on nearest neighbor sites of a hypercubic 
lattice in d space dimensions. The second part describes the interaction of x- 
components of the spins with a magnetic field in transverse (x-) direction. 

Let us first briefly consider the clean model with all Jy = J and all hf = h^. It 
can be viewed, e.g., as a model for the magnetic behavior of LiHoF4 in a transverse 
field [112].+ Without the field, the exchange interaction favors parallel spins pointing 
in positive or negative z-dircction. Thus, for = 0, the ground state is ferromagnetic. 
The transverse magnetic field induces spin flips between the "up" and "down" states 
(25^ = 5+ + 5") and reduces the magnetization in z-direction. For sufficiently strong 
transverse field > /ij!:, this destroys the magnetic long-range order. The quantum 
critical point at was indeed observed experimentally [112]. The partition function 
of the transverse field Ising model (54) in d space dimensions is equivalent to that 
of a classical Ising model in d + 1 space dimensions [12, 113] with imaginary time 
playing the role of the extra dimension. This is an example of the quantum-classical 
mapping discussed in section 2.3. The quantum phase transition of the rf-dimensional 
clean transverse-field Ising model is therefore in the d + 1 dimensional classical Ising 
universality class. 

Quenched disorder can be introduced into the transverse field Ising model by 
making the interaction and/or the transverse field hf a random fimction of position. 
We will only consider the case where all Jij remain positive. (Introducing negative, 
i.e., antiferromagnetic bonds introduces frustration that can lead to qualitatively 
new physics such as spin glass behavior [114,115].) Applying the above-mentioned 
quantum-classical mapping now leads to a disordered classical Ising model. However, 
because the disorder in the quantum model is time-independent, the disorder in the 
equivalent classical model is perfectly correlated in one of the dimensions, viz., the one 
representing the imaginary time direction. Specifically, the one-dimensional random 
transverse- field Ising model (the random transverse Ising chain) maps onto the classical 
McCoy- Wu model [32,33] discussed in section 4.5 when only the interactions Jij are 
random, or on its generalization by Shankar and Murthy [101] when both Jij and hf 
are random. 

Even though the free energy of the random transverse field Ising chain can 
be obtained exactly using transfer matrix methods analogous to those applied to 
the McCoy- Wu model [33, 101], many other properties including the spontaneous 
magnetization and the average spin correlations could not be obtained this way. A 
better understanding of this model was achieved in the early and mid 1990s by two 
methods: (i) a phenomenological optimal fluctuation theory of rare region effects [57] 
and (ii) an asymptotically exact real space renormalization group [34,35]. We will 
discuss these approaches in the next two subsections. 

5.2. Quantum Griffiths singularities: Optimal fluctuation arguments 

In this subsection we develop optimal fluctuation arguments for the rare region effects 
in the random transverse fleld Ising model. These arguments, which were introduced 
by Thill and Huse [57] and applied to both quantum Ising spin glasses [58,59] and the 
random transverse fleld Ising model [39, 116], can be viewed as quantum version of the 

+ The actual interaction in LiHoF4 is of long-range dipolar nature, so the nearest-neighbor interaction 
in the model Hamiltonian (54) constitutes an approximation. 
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arguments given in section 4.5 for the McCoy- Wu model (for an early brief review see 
also reference [117]). 

For definiteness, we consider a situation with random interactions but a 
homogeneous transverse field h^. The probability distribution of the interactions 
is assimied to be binary. P{J) = (1 — p) 5{J — Jq) + pS{J — cJq) with < c < 1. 
Generalizations to other distributions and randomness in the transverse fields are 
straight forward. Consider a rare strongly-coupled spatial region of linear size Lrr 
having only strong bonds. The probability for finding such a region is exponentially 
small in its volume 

w(£rr) ^ exp(-pLRf() . (55) 

where p = — ln(l — p) as before. We are interested in the Griffiths region of the 
disordered phase where the applied transverse field /i^ > but < h^^iean ^^^'^ 
that the rare region is locally ordered. The energy levels of a locally ordered cluster 
can be described as follows [118]. The two lowest states correspond to the states of a 
single effective Ising spin with magnet moment ~ Lf^^i ™ ^.n effective transverse field 
hcs{Ljiji). It can be estimated to be ft.of{(iRR) ^ ft-^ exp[— aLpp] since the "up" and 
"down" states of the effective spin are coupled in perturbation theory of order L"^ in 
h^. The constant a increases continuously with the difference between /i^ clean '^c- 
The energy gap of a rare region of linear size Xrr thus decreases exponentially with 
its volume 

e(LRR) - exp(-aL^R) . (56) 

An analogous result can be obtained by employing the quantum-classical mapping 
and representing the rare region as a one-dimensional classical Ising model [57,59]. 
The exchange interaction in the time-like direction is then proportional to the spatial 
volume of the rare region, Jeff ~ -^rr- Together with the fact that the gap in the 
one-dimensional Ising model depends exponentially on the interaction, this leads to 
(56), see also section 4.5. 

Combining equations (55) and (56) leads to a power-law density of states for the 
low-energy excitations of the rare regions, 

p(e) ~ eP/"-^ = e^l^'-^ (57) 

where the second equality defines the customarily used dynamical exponent z' 
[116]. It continuously varies with disorder concentration or distance from the clean 
critical point. The power-law density of states (57) is obtained here because the 
exponentially small probability for finding a rare region is compensated by the 
exponential dependence of its energy gap on its volume. In contrast, in generic classical 
systems, such as the disordered Ising models of sections 4.1 to 4.3, the gap depends 
on the volume via a power law, and the resulting rare region density of states is 
exponentially small. 

Many results follow from the power-law density of states (57). For instance, a 

region with a local energy gap e has a local spin susceptibility that decays exponentially 
in imaginary time, xioc(''' ^ oo) oc exp(— er). Averaging by means of p yields 

Xr:,(r^^)cxr-'^/^'. (58) 

The temperature dependence of the static average susceptibility is then 

X^c{T)= drxtM^T'^^''-'- (59) 

Jo 
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Figure 6. Left; Scaled probability distribution of the energy gap for a one- 
dimensional random transverse Ising model in the Griffiths region {ho > = 1) 
for different system sizes L. A fit of the low-energy tail to (57) gives z' 1.4 
(from [39]). Right: Specific heat as a function of temperature for the same system. 
A fit to AC oc T''/^' gives 1/z' » 0.75 (from [116]) 



If d < z', the local zero-temperature susceptibility diverges, even though the system 
is globally still in the disordered phase. The contribution of the rare regions to the 
specific heat C can be obtained from 

= Jde p{e) e e-''^ /{I + e"^/^) oc T'^/^'+i (60) 

which gives AC oc T*^/^' . To determine the zero temperature magnetization in a small 

ordering field H in z-dircction wc note that all rare regions with e < H arc (almost) 
fully polarized while the rare regions with e > H have a very small magnetization. 
Thus, 

m.{H) ~ / de p{€) - iJ'^/^ . (61) 
Jo 

Other observables can be determined in a similar fashion. The power-law density of 
states (57) in the Griffiths region of a random transverse field Ising magnet and the 
resulting power-law singularities (58), (59), (60), (61) are called the quantum Griffiths 
singularities or the Griffiths-McCoy singularities. 

The power-law quantum Griffiths singularities predicted by the optimal 
fluctuation arguments above have been verified in numerical simulations. Young and 
Rieger [39] mapped the one-dimensional random transverse field Ising chain onto a 
system of non-interacting fermions which was then solved numerically. An example 
of the resulting density of states in the Griffiths region is shown in the left panel of 
figure 6. The low-energy tail of this distribution shows power-law behavior. A fit to 
the predicted form (57) gives 2;' w 1.4 for the parameter set studied. The right panel 
of figure 6 shows the temperature dependence of the specific heat for the same system 
as calculated by Young [116]. As predicted after (60), this dependence is of power-law 
type, and a fit gives the exponent z' « 4/3 in good agreement with the value from 
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the density of states. These numerical simulations have later been extended to surface 

quantities and to the ferromagnetic phase [119] as well as to nonlinear susceptibility, 
higher excitations, and the energy autocorrelation function [120]. 

5.3. Strong- disorder renormalization group and infinite-randomness critical point in 
one dimension 

While the optimal fluctuation arguments presented in the last subsection provide a 
simple qualitative understanding of the behavior of random transverse field Ising model 
in the Griffiths region, they do not give quantitative results and they also cannot be 
applied to the quantum critical point itself. The behavior of the random transverse 
field Ising model was therefore not fully understood before Fisher [34, 35] applied a 
novel real-space or rather energy-space renormalization group method to the problem. 
This method which is now often called the "strong-disorder renormalization group" 
has recently been reviewed in detail by Igloi and Monthus [121]. Therefore we will 
only briefly explain its basic idea and then discuss the results. 

The strong-disorder renormalization group was introduced by Ma, Dasgupta and 
Hu [36] for the random antiferromagnetic spin chain. Since then it has been used 
extensively for various disorder problems including random quantum spin chains 
[38, 122-133] and ladders [134-137] as well as random walks [138-140] and other 
nonequilibrium systems [41, 42, 141, 142]. The main idea of the strong-disorder 
renormalization group is to take the strongest coupling in the system and find the 
ground state of the corresponding part of the Hamiltonian exactly. The coupling to 
the rest of the system is treated perturbatively. One then neglects the excited states 
involving the strong coupling and derives a new effective Hamiltonian with the number 
of degrees of freedom reduced. This basic renormalization group step is repeated ad 
infimum. 

We now follow reference [35] and apply this idea to the random transverse 
field Ising model (54) in one dimension. We choose first the largest of the set of 
couplings, ft = max(ft,f , Jij). If the largest couphng is a field, say /if, the unperturbed 
Hamiltonian is — /ifSf with ground state | ^2) and excited state | ^2) separated by 
a gap /if. The coupling of 52 to the rest of the system — Ji2 5'f 5*1 — J23<5'|S'| is treated 
in second order perturbation theory. This yields an effective interaction of the form 
—J-isSfS^ with the effective exchange constant 

Ji3 « J12 J23//if • (62) 

We now throw out 5*2 completely by neglecting its excited state. This loads to a new 
spin chain with one less spin and a new coupling J13 < fJ. If the largest coupling is 
a bond, say, J23, the corresponding unperturbed Hamiltonian is — J23S'|5'| with two 
degenerate ground states | ft) and | J, J,). The transverse fields /if and /if can now be 
treated perturbatively yielding an effective field 

/if « hyj23 (63) 

inducing coherent flips of the cluster 5*2-1-53. We now throw away the exited 
(antiparallel) states of this cluster and treat the cluster as an effective spin whose 
moment is the sum of the moments of 5*2 and 53. Again, we have a new spin chain 
with one less spin degree of freedom and all couplings smaller than fl. 

The repeated application of this basic renormalization group step constitutes a 
aggregation and annihilation process of spin clusters. When the strongest coupling is 
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a field, the corresponding cluster is annihilated, and when it is a bond, the clusters 
that it connects arc aggregated into one cluster. At each renormalization group step, 
the clusters represent sets of original spins that at correlated at the current cutofi: 
energy O. In the paramagnetic phase, annihilation dominates for O ^ 0, and no 
large clusters arc generated. In the ferromagnetic phase, the aggregation dominates 
for n ^ generating arbitrarily large clusters until, at f2 = an infinite cluster will 
have formed. The quantum critical point is the point where such an infinite cluster 
first appears [26, 126]. The multiplicative structure of the recursion relations (62) and 
(63) is very important because it establishes an exponential relationship between the 
length scale and the energy scale: In each RG step, the lengths of clusters or bonds 
arc added while their couplings are multiplied to determine the effective renormalized 
values. Thus, the length scale L of the clusters and efi'ective bonds scales with the 
logarithm of the energy scale 

L'^ ~ ln(fio/f2) (64) 

where fio is a microscopic energy scale and the critical exponent tp turns out to be 
ip = 1/2. Analogously, the typical magnetic moment of a cluster increases as 

/X - (ln(r!o/0))^ (65) 

with the exponent cj) = (1 + •\/5)/2 equal to the golden mean. Thus, the 
random transverse field Ising model displays activated dynamical scaling rather than 
conventional power-law scaling. 

Using this approach. Fisher [34,35], derived renormalization group equations for 
the probability distributions of the logarithms of the couplings. In J and In/i^. The 
resulting renormalization group flow displays very interesting and unusual properties. 
At the critical point which occurs at [In J]av = [ln/i,^]av where [. . .]av is the disorder 
average [101], the probability distributions of In J and \nh^ broaden without limit with 
decreasing energy scale D,. This broadening is crucial for the renormalization group 
to work: In the limit of infinitely broad distributions, the recursion relations (62) and 
(63) become asymptotically exact for ^ 0, since in each renormalization group step 
the ratio between the dominating coupling and those treatc^d perturbatively diverges. 
Because of the diverging widths of the probability distributions of In J and \nh^, the 
critical fixed point in the random transverse field Ising model is also called an infinite 
randomness fixed point. 

Off criticality, the distance from the critical point can be conveniently defined as 

^ ^ [In/llav- [lnJ]av .ggs 

var(ln/i^) + var(ln J) ^ ' 

where var(.'j;) is the variance of the probability distribution of x. The so defined r is 
invariant under renormalization. In the disordered (paramagnetic) phase, there will 
be a length scale ^ at which the renormalization flow qualitatively changes. On scales 
larger than ^ almost all effective fields are larger than the effective couplings. Upon 
further renormalization, only fields will decimated, leading to longer and weaker bonds 
but not to longer clusters. Thus, ^ is the correlation lengths. It is found to scale like 

^ ~ r"'' (67) 

where the exponent v = 2 saturates the Harris criterion du > 2 [25]. Analogously, 
in the ordered (ferromagnetic) phase for r small and negative, there will be a length 
scale ^ [r]"" beyond which (almost) all interactions are larger than the fields. Upon 
further renormalization, all clusters will be joined together to form the infinite cluster 
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at f2 = 0. Thus, excitations from the ferromagnetic state exist only up to the length 

scale ^ which is thus the (connected) correlation length in the ordered phase. The three 
exponents 0, tp and p completely determine the properties of the infinite-randomness 
critical point. 

The imusual rcnormalization group flows lead to correspondingly unusual scaling 
behavior of observables at the quantum phase transition. At zero temperature, the 
magnetization in a small ordering magnetic field H in ^;-direction can be determined 
by stopping the rcnormalization at an energy scale ^Ih = H and analyzing the cluster 
distribution at that scale.* This leads to the scaling form 

m (r, \n{Ho/H)) = b^^~'^ m [rb^/-" MHo/H)b-^) (68) 

where b is an arbitrary scale factor, and Hq is a nonunivcrsal constant. This equation 
again reflects the activated nature of the dynamical scaling; the energy scale set by 
H enters the scaling relations in logarithmic form rather than in power-law form. At 
criticality, r = 0, we can set b ^ [\n{Ho/H)]^/'l' and obtain 

m.{H) - [IniHo/H)]'^'-'^^'^ (69) 

The exponent evaluates to (p — d/ip = </> — 2 = (v^— 3)/2. Setting H = and b = r~'' 
we can obtain the spontaneous magnetization in the ordered phase (r < 0) 

m(r) ~ \rf (70) 

with /3 = {d — (j)'il})u = 2 — (j)= (3 — \/5)/2. Temperature dependent quantities can be 
determined in a similar fashion. At criticality and zero ordering field -ff , the magnetic 
susceptibility behaves as 

x-^Mnrnff-'^'^ (71) 

where Tq is a nonunivcrsal constant, and the specific heat is given by 

~ [ln(To/r)]-(i+'^/^) . (72) 

The infinite-randomness character of the critical point also has interesting 
consequences for the behavior of the spin-spin correlation function C{x) = 
{SfSf_^_^). Because the coupling constant distributions broaden without limit under 
rcnormalization, C{x) is also very broadly distributed. As a result typical and average 
correlations behave very diff'erently. At the critical point, typical spin pairs are never 
in the same rcnormalization cluster and thus have only weak correlations that fall off 
as a stretched exponential 

\nCtyp{x) - -x"^ . (73) 

In contrast, the average correlation function is dominated by rare spin pairs on the 
same cluster that have atypically large correlations of order unity. It is thus much 
larger than the typical one and falls off as a power law 

C^^ix) ^ x-^^'^-^^K (74) 

This discrepancy between typical and average properties is one of the main 
characteristics of an infinite-randomness critical point. Fisher and Young [143] 
numerically studied the end-to-end correlations in an open one-dimensional random 

* More precisely, the field energy scale is Hfi where n is the typical cluster moment at scale CIh- 
This introduces an extra logarithmic correction. 



Rare region effects at classical, quantum, and nonequilibrium phase transitions 35 

transverse-field Ising model and found good agreement with the prediction of the 
strong-disorder rcnormalization group. 

The strong-disorder renormalization group can also be used to analyze 
thermodynamic observables and correlation functions off criticality. In the Griffiths 
region of the disordered phase. Fisher [35] found strong power-law quantum Griffiths 
singularities of the form predicted by the optimal fluctuation arguments presented in 
the last subsection. Moreover, he found an explicit expression for the nonuniversal 
dynamical exponent z' . Close to the true quantum critical point, i.e. for r <C 1 it is 
given by z' = 1/ 1 2r | . The discrepancy between typical and average correlations persists 
in the Griffiths region. Typical correlations fall off exponentially with a correlation 
length ^typ that diverges as ^typ ^ [rj^'^typ with i^typ = 1- In contrast, the average 
correlations decay with the much larger true correlation length ^ ~ with u = 2. 
An analogous Griffiths region with power-law singularities was also identified in the 
ordered phase. 

Igloi et al. [144, 145] extended Fishers analytical solution into the off-critical 
region and obtained asymptotically exact results for dynamical quantities including 
the dynamical exponent z' . Static quantities or spatial correlations, however, cannot 
be obtained exactly off criticality because the coupling constant distributions do not 
broaden without limit under renormalization. Rieger and Igloi [146] studied the 
influence of power-law spatial correlations of the random couplings on the transition. 
They found the infinite-randomness scenario still valid, but sufficiently long-ranged 
correlations yield new exponents that enhance both the critical and the quantum 
Griffiths singularities. 

5.4- Higher dimensions and other generalizations 

The discovery of an infinite-randomness critical point in the one-dimensional random 
transvcrse-flcld Ising model lead to a surge of activity in this area starting in the 
mid 1990's. Early on, an important question was whether or not the infinite- 
randomness critical point is special to one dimension, or whether it also occurs in 
higher-dimensional systems. Senthil and Sachdev [118] considered a bond-diluted 
transverse field Ising model in d > 2. A sketch of the phase diagram of a general 
diluted quantum magnet is shown in figure 7. The diluted transverse field Ising model 
has two quantum phase transitions, separated by a multicritical point at the lattice 
percolation threshold Pc and a transverse field /i^ = h* . The transition at dilutions 
p < Pc (transition "a" in figure 7) is expected to be in the universality class of the 
generic random transverse-field Ising model. Senthil and Sachdev focused on the 
transition "b" across the percolation threshold of the lattice. Using the well known 
properties of critical percolation clusters [147], they showed that this quantum phase 
transition shares many properties with the infinite-randomness critical point in the 
one-dimensional random transverse Ising model. Specifically, the dynamical scaling 
is activated and the magnetization in a small ordering field obeys the scaling form 
(68) with ip = Df, <p = 1 and v = Vc where Df is the fractal dimension of the 
critical percolation cluster and Vc is the correlation length exponent of the lattice 
percolation problem. The transition is accompanied by strong power-law quantum 
Griffiths singularities both above and below pc- All these predictions were verified by 
quantum Monte-Carlo results of Ikegami et al. [148]. 

To study the generic random transverse field Ising transition in dimensions d>2, 
Motrunich et al. [26] applied a strong-disorder renormalization group to the two- 
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Figure 7. Schematic phase diagram of a general diluted quantum magnet with 
impurity concentration p, temperature T, and quantum fluctuation strength g 
(for the transverse- field Ising model g ^ h'^ . There is a quantum multicritical 
point (big dot) at (pci9*)- 

dimensional random transverse-field Ising model. In dimensions larger than one, 
the renormalization group calculation cannot be carried out analytically because the 
renormalization group does not preserve the lattice structure but generates longer 
and longer ranged couplings. Therefore, Motrunich et al. [26], implemented the 
renormalization procedure numerically. For two dimensions, they found that the 
renormalization group flow for stronger disorder is indeed towards even stronger 
disorder, i.e., towards an infinite-randomness fixed point, as in one dimension. As 
a result, the same scaling scenario as outlined in section 5.3 also applies to the two- 
dimensional case, but with different exponent values. Motrunich et al. estimated 
"ip K, 0.42, (j) « 2.5, and v ~ 1.1. Other numerical implementations [149, 150] of 
the strong disorder renormalization group found similar values. Preliminary results 
for three dimensions point to the same scenario, but reliable estimates for the 
exponents have not yet been obtained. The results of the numerically implemented 
strong-disorder renormalization group are in very good agreement with Monte-Carlo 
simulations of a two-dimensional random transverse field Ising model [40, 151]. 

The infinite-randomness critical fixed points found in the random transverse field 
Ising model control more than just Ising ferromagnetic quantum phase transitions. As 
was first pointed out by Senthil and Majumdar [152], Potts models or any random 
quantum systems with continuous quantum phase transitions at which a discrete 
symmetry of a non-conserved order parameter is broken, will have the same critical 
behavior as the Ising case. This applies even to systems that are frustrated on 
microscopic length scales, such as quantum Ising spin glasses. Because the coupling 
constant distributions become infinitely broad at the infinite-randomness fixed point, 
in any frustrated loop the weakest interaction will be infinitely weaker than the others 
and can thus be neglected. Thus, for Ising symmetry, frustration becomes irrelevant 
at the infinite-randomness critical point. 
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5.5. Continuous symmetry quantum magnets 

In the last subsections we have seen that the qualitative zero-temperature behavior of 
random quantum Ising systems docs not depend much on dimensionality: There is an 
infinite-randomness critical point separating a conventional paramagnetic phase from 
a phase with long-range order (of course, the critical exponent values at this critical 
point do depend on dimensionality). 

For continuous- symmetry quantum spin systems such as quantum Heisenberg or 
XY magnets, the situation is drastically different. In one space dimension, these 
systems do not show magnetic; long-range order. Instead they can be in a variety 
of exotic phases that will be briefiy discussed in section 5.7. In staying with the 
main topic of this review, the present subsection is focussed on transitions between 
a disordered and a long-range ordered phase that can occur in continuous symmetry 
quantum magnets for d> 2. 

Let us start by considering a two-dimensional 5=1/2 Heisenberg quantum 
antiferroniagnet with random site dilution. This problem is of direct relevance for 
antiferromagnetic layered cuprates with nonmagnetic impurities [153-156]. In the 
clean limit, the ground state of this system displays magnetic long-range order. 
However, the zero-point quantum fluctuations reduce the staggered magnetization by 
about 40% compared to the classical value (see, e.g., [157]). Upon randomly removing 
sites (or bonds) with probability p, the tendency towards magnetic order decreases, 
but the location and character of the phase transition were initially controversial. 
Early theoretical investigations [158-161] suggested that magnetic long-range order 
is destroyed before the impurity concentration p reaches the classical (geometric) 
percolation threshold pc ~ 0.407 [147] of the lattice. Kato et al. [162] found evidence 
of the critical dilution coinciding with percolation threshold Pc but argued that the 
infinite percolation cluster is quantum critical, with power-law spin correlations. 
However, using extensive quantum Monte-Carlo simulations, Sandvik [163,164] showed 
that the infinite percolating cluster at Pc displays long-range order. The transition 
at Pc is thus of percolation type. If the quantum fluctuations are enhanced to 
induce a quantum phase transition at impurity concentrations p < pc, the zero- 
temperature phase diagram takes the general form shown in figure 7. However, at finite 
temperatures, no long-range order exists because of the Mermin- Wagner theorem [165]. 

An obvious and experimentally relevant way to increase the quantum fluctuations 
is via an interlayer coupling in a bilayer Heisenberg quantum antiferromagnet as 
depicted in the inset of figure 8. The spins in each 2d layer interact via nearest 
neighbor exchange J||, and the interplane coupling is J± with the ratio g ^ J±/J\\ 
controlling the strength of the quantum fluctuations. The clean version of this model 
has been studied extensively [166-169]. For J± » Jy, neighboring spins from the two 
layers form singlets, and the ground state is paramagnetic. In contrast, for Jy ^ J± 
the system develops Neel order. Both phases are separated by a quantum phase 
transition at J±/J\\ ~ 2.525. 

Quenched disorder can be introduced by randomly removing spins. We emphasize 
that in order to achieve random-mass type disorder, pairs (dimers) of adjacent spins, 
one from each layer have to be removed. In this way, no random Berry phases 
arise, because the Berry phase contributions from the two spins of each unit cell 
exactly cancel [12,170]. In contrast, for site dilution, the physics would completely 
change: The random Berry phases (which have no classical analogue) are equivalent 
to impurity-induced moments [171], and those become weakly coupled via bulk 
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Figure 8. Phase diagram [68, 69] of tiie diluted bilayer Heisenberg 
antiforromagnet, as function of J^/Jy and dilution p. The dashed line is the 
percolation threshold, the open dot is the multicritical point. Inset: The model: 
Quantum spins (arrows) reside on the two parallel square lattices. The spins in 
each plane interact with the coupling strength Jy. Interplane coupling is J±. 
Dilution is done by removing dimers (from [70]). 



excitations. For all p < Pp the ground state of the system shows long-range order, 
independent of J±/J^^l Thus, site dilution actually changes the bulk phases, i.e., it 
introduces disorder effects more complicated than random-mass disorder. These are 
beyond the scope of this review. 

The Hamiltonian of the dimer-diluted bilayer Heisenberg model is: 



H = Ji 



E 

a=l,2 



1,1 



(75) 



and 6^=0 {ei=l) with probability p (1 — p). The zero-temperature phase diagram of 
this model has been studied in detail by Sandvik [68] and Vajk and Greven [69]; it is 
shown in Fig. 8. For small J±, magnetic order survives up to the percolation threshold 
Pp ~ 0.4072, and a multicritical point exists at p = pp and J±/ J\\ ~ 0.16. Thus, the 
dimer-diluted bilayer Heisenberg antiferromagnet has two quantum phase transitions, 
the generic transition for p < pc (corresponding to transition "a" in the generic phase 
diagram, figure 7) and a quantum percolation transition at p = Pc,Ji_ < 0.16J|| 
(corresponding to transition "b" in figure 7) . 

The low-energy properties of bilayer quantum antiferromagnets are represented 
by a (2-|-l)-dimensional 0(3) quantum rotor model with the rotor coordinate fii 
corresponding to Si.i — Si.2 and the angular momentum Li representing Si,i -|- Si, 2 
[12, 170]. This quantum rotor model in turn maps onto a 3d classical Heisenberg model 
in which the disorder is perfectly correlated in the imaginary time direction. At zero 
temperature, the rare regions are thus equivalent to classical Heisenberg magnets that 
are infinite in the time-like direction but finite in the two space directions. 

A first qualitative understanding of the rare region effects can again be obtained 
from optimal fluctuation arguments similar to that used in section 5.2. The probability 
for finding a rare region of linear size Lrr, is exponentially small in its volume, 
wiLji^) ^ exp(— pLpfj) with d — 2. However, the energy gap of a locally ordered 
rare region depends on its spatial volume via a power law (rather than exponentially 
as in the Ising case) as can be obtained from a renormalization group analysis of the 
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(0+1) dimensional Heisenberg model at its strong coupling fixed point or by an explicit 
large- iV calculation [67]. These arguments thus predict the rare region effects in the 
Heisenberg case to be exponentially weak and the critical point to be conventional 
with power-law dynamical scaling. 

Sknepnek, Vojta and Vojta [70] performed large scale Monte-Carlo simulations 
of a classical three-dimensional Heisenberg model with linear defects to test these 
predictions. They focused on the generic transition for p < Pc and analyzed the 
scaling behavior of the so-called Binder ratio of the order parameter moments, given 
by 

(Iml-*) 1 , , 

\l I / ^7g^ 



1 



3(|m|2)2 

where [. . .]^^ denotes the disorder average and (. . .) is the Monte-Carlo average for 
each sample. This quantity has scale dimension 0. Thus, its finite-size scaling form is 
given by 

gav=gcirL'/-',Lr/L') or (77) 
5a. =5A(rLl/^log(L0/L'^) (78) 
for conventional scaling or for activated scaling, respectively. Here L and Lt are the 
linear system sizes in the space and time-like directions. Two important characteristics 
immediately follow from these scaling forms: (i) For fixed L, g^v has a peak as a 
function of L^- . The peak position i™ax j^iarks the optimal sample shape, where the 
ratio Lr/L roughly behaves like the corresponding ratio of the correlation lengths in 
time and space directions, ^r/C- -^t the critical temperature Tc, the peak value g"^^^ 
is independent of L. Thus, for power law scaling, plotting gav vs. Lr/L^°'^ at Tc 
should collapse the data, without the need for a value of z. In contrast, for activated 
scaling the gav data should collapse when plotted as a fimction of log(LT-)/ log(L™^'^). 
(ii) For samples of the optimal shape {Lt = L™"^), plots of gav vs. temperature for 
different L cross at T^. 

A typical result of these calculations is presented in figure 9. This figure shows the 
Binder ratio at the critical point for a system with impurity concentration p = 1/5 
and various system sizes L and Lr- As can be seen in the main panel, the data 
scale very well when analyzed according to conventional power-law scaling while 
the inset shows that they do not scale at all when plotted according to activated 
scaling. Analogous data where obtained for impurity concentrations 1/8, 2/7, and 
1/3. Sknepnek et al. [70] thus concluded that the quantum critical point in the dimer- 
diluted bilayer Heisenberg antiferromagnet exhibits conventional power-law scaling 
with exponentially weak rare region effects, in agreement with the optimal fluctuation 
arguments above. They also determined the critical exponents at this quantum 
critical point and found universal (i.e., disorder-independent) values oi z = 1.310(6), 
V = 1.16(3), P/iy = 0.56(5), and j/u = 2.15(10). 

Sandvik [68] and Vajk and Greven [69] studied the scaling behavior at the 
multicritical point a.t p = Pp and J±/J\\ « 0.16. They found that this point also 
exhibits conventional power-law scaling. Their value of the dynamical exponent 
2 « 1.3 agrees within the error bars with that of the generic transition discussed 
above. Whether or not the dynamic exponents at the generic transition and at the 
multicritical point are indeed identical and whether or not this agreement is accidental 
is not resolved so far. 

More recently, Vojta and Schmalian [172] considered the quantum phase 
transition across the percolation threshold (transition "b" in the generic phase 
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Figure 9. Upper panel: Binder ratio gav as a function of Lt for various L 
(p = i). Lower panel: Power-law scaling plot gav/gl^^ vs. Lt/L™"^ Inset: 
Activated scaling plot gav/gZ?"^ vs. y = \og(L-r) / \og(L^'^'') (from [70]). 



diagram, figure 7). Combining the well known properties of the lattice percolation 
transition [147] with a scaling analysis of finite-size quantum Heisenberg clusters at the 
stable low-energy fixed point, they derived a complete scaling theory of this transition. 
In contrast to the corresponding quantum Ising case [118], the critical behavior is of 
power-law type. The dynamical exponent is identical to fractal dimension of the 
percolation cluster, z = Df, and the correlation length exponent is identical to that 
of the lattice percolation problem, = In fact, these results apply not only to the 
dimer-diluted bilayer model but general diluted Heisenberg quantum antiferromagnets. 
The conclusion z = Df has been confirmed by Monte Carlo simulations of single- 
layer [173] and bilayer [68] models. 

We can thus conclude that all quantum phase transitions of the dimer-diluted 
bilayer quantum Heisenberg antiferromagnet (the generic transition for p < pc, the 
percolation transition at pc and the multicritical point) exhibit conventional power- 
law scaling and exponentially weak rare region effects in agreement with the simple 
extremal statistics arguments. This conclusion is also supported by preliminary results 
of a strong-disorder renormalization group calculation for the Heisenberg case [26] that 
show that the infinite-randomness fixed point is unstable and the renormalization 
group flow is towards smaller disorder strength. Lin et al. [174] used the strong- 
disorder renormalization group to study the stable low-energy (as opposed to critical) 
fixed points of random Heisenberg models in two and three dimensions. They found 
conventional power-law scaling and no infinite randomness behavior. This result has 
recently been confirmed by quantum Monte-Carlo simulations and exact numerical 
diagonalizations [175]. 
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5.6. Itinerant quantum magnets 

In this subsection, we turn from systems of localized spins to quantum phase 
transitions in itinerant (metallic) spin systems. This field was pioneered by Hertz 
[176]. Starting from a microscopic theory of interacting electrons, he derived 
Landau-Ginzburg- Wilson order parameter field theories for the ferromagnetic and 
antiferromagnetic quantum phase transitions of both clean and dirty metals. He 
then studied these Landau-Ginzburg- Wilson theories using renormalization group 
methods. Later, Millis [177] worked out the experimentally important finite- 
tempc;ratur(^ behavior in the vicinity of the quantum critical point. In tlw^ c;ase of 
the ferromagnetic transition, the specific conclusions of the Hertz-Millis theory turned 
out to be incorrect because generic scale invariance in the of form additional fermionic 
soft modes invalidates the Landau-Ginzburg- Wilson approach [14,18 20]. This was 
briefly discussed at the end of section 2.3. The effective long-range interactions created 
by the additional soft modes will likely modify the rare region effects, too. However, 
these extra complications at the dirty ferromagnetic transition are presently not fully 
understood, and therefore we will not discuss them in detail. 

Physically, a crucial difference between itinerant magnets and systems of localized 
spins is that magnetic excitations are damped in the former while they are undamped 
in the latter. This is caused by the coupling between the magnetic modes and 
the gapless particle-hole excitations in the metal. Technically, it is reflected in the 
frequency dependence of the Gaussian term of the Landau-Ginzburg- Wilson theory. 
For definiteness, let us consider the itinerant antiferromagnetic transition. The 
Landau-Ginzburg- Wilson free energy functional of the clean transition reads [176]| 



Here, r(x, r, y,T') is the bare two-point vertex whose Fourier transformation is 



with r being the (bare) distance from criticality and cOn bosonic Matsubara frequencies. 
We have suppressed constant prefactors. The nonanalytic frequency dependence 
proportional to \LOn\ in the dynamical part of T reflects the overdamping of the 
dynamics due to the coupling of the order parameter to fermionic particle-hole 
excitations. In contrast, undamped dynamics would lead to an term (see also 
equation (11)). 

Weak, random-mass disorder can be introduced by making r a random function 
of position, r — > r -|- 5r{x) [176, 179]. The rare regions in this system are large spatial 
regions where the local r is smaller than its average value. Let us consider a single 
such region that is locally in the ordered phase, while the bulk is still in the disordered 
phase. The order parameter on such a region is correlated in space, i.e., it behaves 
coherently over the whole island. The crucial difference between itinerant magnets 
and the localized spin systems considered in previous sections is in the dynamics of 
the rare region. Without damping, the order parameter slowly fluctuates giving rise 
to the quantum Griffiths effects discussed in section 5.2. Damping or dissipation will 
potentially hinder the dynamics of the rare region. To study this phenomenon, Millis, 
Morr and Schmalian [61,62] explicitly calculated the tunneling rate of a locally ordered 
rare region in an itinerant Ising magnet. They found that for sufficiently large rare 

I This simple bosonic Landau-Ginzburg Wilson theory only applies for dimensions d > 2. In two 
dimensions, fermionic soft modes lead to singular vertices even for the antiferromagnet [178]. 




(79) 
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regions, the tunneling rate vanishes. This means, these rare regions completely stop 
to tunnel, and their order parameter becomes static. The same result can also be 
obtained from quantum- to-classical mapping [60]. In the equivalent classical system, 
a rare region corresponds to a quasi-one-dimensional Ising model which is of finite size 
in the space directions but infinite in the time-like direction. The linear frequency 
dependence in the two-point vertex F is equivalent to a long-range interaction in 
imaginary time of the form (r — t')~^. Each rare region is thus equivalent to an one- 
dimensional Ising model with a 1/r^ interaction. This model is known to have a phase 
transition [180,181]. Thus, true static order can develop on those rare regions which 
are locally in the ordered phase. § Very recently, Schehr and Rieger [182] analyzed 
a related problem, viz. a random transverse field Ising model coupled to an Ohmic 
dissipative bath, by a version of the strong-disorder renormalization group. They 
arrive at an analogous conclusion: Sufficiently large, locally ordered clusters freeze, 
and develop a static order parameter. 

As a result, the phenomenology of rare region effects in an itinerant Ising magnet 
is very different from that of a localized magnet. Because the large rare regions do 
not fluctuate the low-energy behavior docs not display power-law quantum Griffiths 
effects, contrary to earlier predictions that did not correctly take into account the 
dissipation [183]. Instead, the global phase transition is smeared [60] by the same 
mechanism as the transition in a classical Ising model with planar defects discussed 
in section 4.6: Once static order has developed on a few isolated rare regions, an 
infinitesimally small interaction or an infinitesimally small symmetry-breaking field 
arc sufficient to align them. Consequently, a macroscopic order parameter develops 
inhomogeneously, with different spatial part of the system ordering at different values 
of the control parameter. 

The leading zero-temperature behavior in the tail of the smeared quantum phase 
transition can be determined by optimal fluctuation arguments similar to section 4.7. 
Here, we just summarize the results for dilution- type (Poisson) disorder, i.e., Sr{x) = 
everywhere except on randomly distributed finite-size islands (impurities) of spatial 
density p where (5r(x) = W > 0. More details can be found in reference [60]. The 
probability w for finding a region of linear size Lrr devoid of any impurities is given by 
w ^ exp(— pipj^) (up to pre-exponential factors). Such a rare region develops static 
order at a distance rc(iRR) < from the clean critical point. Finite size scaling yields 
|''c(iRR)| ~ where is the finite-size scaling shift exponent of the clean system. 
Thus, the probability for finding a rare rcigion which becomes critical at Vc is given 
by w{rc) ~ exp(— i? |rc|~'^/'^) for r 0—. The total order parameter m is obtained 
by integrating over all rare regions which are ordered at r, i.e., all rare regions having 
tc > t. This leads to an exponential tail of m as a function of the distance r from the 
clean critical point: 

m{r) ^ oxp{-B\r\-'^/'t') (for r ^ 0-) . (81) 

This is identical to the classical result (47) if one replaces the number of uncorrelated 
dimensions rfran by the space dimensionality d. 

We now turn to the experimentally important finite-temperature behavior in the 
vicinity of the smeared quantum phase transition. A schematic phase diagram is 
shown in figure 10. At any finite temperature, the static order on the rare regions 

§ At the first glance, this result seems to contradict basic statistical mechanics that tells us that 
spontaneous symmetry breaking cannot occur in a finite-size system. The resolution lies in the fact 
that the damping is caused by the underlying electronic system which must be infinite to be gapless. 




Figure 10. Schematic phase diagram in the vicinity of a smeared quantum phase 
transition. For details see text. 



is destroyed, and a finite interaction of the order of the temperature is necessary 
to align them. This means a sharp phase transition is recovered. To estimate the 
transition temperature we note that the interaction between two rare regions depends 
exponentially on their (typical) spatial distance |x| which itself depends exponentially 
on r. This leads to a double-exponential dependence of the critical temperature on r 

log(-alogTe)^ |r|-^/^ . (82) 

In other words, the long-range ordered phase develops an double-exponential tail 
towards to clean quantum critical point at r — 0. (Note that the inclusion of a 
long-range spatial interaction such as the RKKY interaction would transform this 
into a simple stretched exponential. Moreover, weak competing interactions not 
contained in the model (79) may induce spin glass order in the tail region.) Above 
this very low temperature different rare regions are not aligned anymore but act as 
independent classical moments leading to super-paramagnetic behavior with a Curie- 
type susceptibility [62]. 

What happens at even higher temperatures depends on nonuniversal microscopic 
parameters and is subject of an ongoing debate. If the damping of the magnetic modes 
by the electrons is weak, i.e., if the prefactor of the |a;„|-term in the Gaussian vertex 
(80) is parametrically small compared to other scales, a subleading term (which 
describes undamped dynamics) becomes dominating at higher energies. This gives rise 
to a crossover temperature T* . Power-law quantum Griffiths behavior can occur in a 
transient temperature window above T* but below a microscopic cutoff scale (like the 
Fermi temperature). However, if the damping is not weak, T* should roughly coincide 
with the microscopic cutoff, and quantum Griffiths effects should not be observable. 
For more details see the discussion on experiments in section 7.2. 

So far, our exploration of itinerant quantum magnets has focused on the case 
of Ising symmetry where we found a smeared phase transition and static (frozen) 
rare regions. For continuous spin symmetry, the results are markedly different, as 
can be readily understood from quantum-to-classical mapping arguments [67] . A rare 
region in an itinerant magnet with continuous spin symmetry is equivalent to a quasi- 
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one-dimensional classical 0{N) model {N > 2). The interaction in the time-like 
direction has the form (t — r')"^ as in the Ising case. Onc-dimcnsional continuous- 
symmetry 0{N) models with interaction do not have a phase transition, but 
they are exactly at their lower critical dimension [184-186]. Therefore, in the 
quantum problem, an isolated rare region of linear size inn cannot independently 
undergo a phase transition. Its energy gap depends exponentially on its volume, 
({Lkr) ~ exp(— aL^j^). Equivalently, the susceptibility of such a region diverges 
exponentially with its volume. Combining this exponential law for the gap with 
an the exponentially small probability for finding a rare region leads to exactly the 
same power-law quantum Griffiths singularities as in the random transverse field Ising 
model [67]. Loh ct al. [187] studied a related problem, viz., a single spin-S impurity 
with XY rotational symmetry in a metal. Their results are completely compatible 
with the quantum Griffiths scenario outlined above. 

This scenario is correct as long as different rare regions are effectively independent. 
Dobrosavljevic and Miranda [73] pointed out that in realistic metallic systems the 
long-range RKKY interaction leads to an additional damping of the rare regions that 
destabilizes the quantum Griffiths behavior and leads to static frozen rare regions even 
in the Heisenberg case. However, in contrast to the Ising case, the freezing requires an 
interaction between the rare regions and is therefore expected to occur at extremely 
low temperatures only. 

5. 7. Quantum Griffiths effects in other systems 

In accordance with the main theme of this review, our discussion of quantum Griffiths 
effects so far has been restricted to order-disorder quantum phase transitions between 
conventional phases, i.e., we have assumed that the phases themselves are not 
qualitatively changed by the quenched disorder. In the present subsection we will 
briefly mention, without any pretense of completeness, other examples of quantum 
rare region effects. 

One-dimensional Heisenberg quantum antiferromagnets do not display long-range 
order even at zero temperature. The behavior of the clean chain depends on the spin 
value [188, 189]: For integer spin, there is a Haldane gap for magnetic excitations and 
short-range magnetic order. In the case of half-integer spin, there is quasi long-range 
spin order with algebraically decaying correlations and no gap. This suggests that the 
half-integer chains are more susceptible towards weak disorder. Indeed, the ground 
state of the random antiferromagnetic spin-1/2 chain is a random-singlet state, that is 
controlled by an infinite-disorder fixed point [38] . Numerical simulations mostly agree 
with this prediction [190, 191]. In contrast, the dimerized spin-1/2 chain has a gap in 
the absence of quenched disorder. As a result, the dimerized chain is stable against 
disorder [123] and the spin correlations remain short-ranged (random dimer phase). 
For strong disorder or weak dimerization, the random dimer phase is gapless, with a 
susceptibility that diverges as some nonuniversal inverse power of temperature. This 
phase can be understood as a Griffiths phase belonging to the critical point at zero 
dimerization. 

The clean antiferromagnetic spin-1 chain has a Haldane gap. Consequently, 
it is stable against weak disorder. Strong-disorder renormalization group studies 
[124, 126, 127, 129] predict a random-singlet phase for sufficiently strong disorder. It 
is separated from the gaped phase by a gapless region of Griffiths singularities, the 
so-called gapless Haldane phase. While the numerical simulations [192-196] support 
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the overall structure of the phase diagram, they give somewhat conflicting results on 
what phases can be reached for which bare disorder distribution and on the location 
of the critical points. For the spin-3/2 chain, even the question of the phase diagram 
has not yet been settled. A strong disorder renormalization group study by Refael 
et al. [130] gave two random singlet phases with different effective spins {Sdi = 3/2 
for strong disorder and SeS = 1/2 for weak disorder). A different implementation 
of the renormalization group by Saguia et al. [131] gives a random singlet phase 
only for strong disorder. Weak disorder is an irrelevant perturbation. The latter 
result also agrees with recent numerical simulations [197]. In addition to spin chains, 
disordered spin ladder systems also display random singlet, Griffiths and other exotic 
phases [132,134-137]. 

Interesting quantum Griffiths effects have also been identified in inhomogeneously 
diluted quantum antiferromagnets. As was discussed in section 5.5, in a 
(homogeneously) diluted two-dimensional Heisenberg quantum antiferromagnct, the 
ordered phase persists to the classical percolation threshold of the lattice, and 
the quantum phase transition that eventually destroys the long-range order is of 
percolation type. Yu et al. [173,198] introduced inhomogeneous bond dilution by using 
different occupation probabilities for different types of bonds ("dimer" and "ladder" 
bonds). In this way, they increase the quantum fluctuations and the local tendency 
towards valence-bond solid order. If the two occupation probabilities are sufficiently 
different, magnetic long-range order is destroyed before the classical percolation 
threshold is reached. In the resulting quantum disordered phase, the low-temperature 
susceptibility diverges algebraically with a non-universal exponent. This phase can be 
understood as a quantum Griffiths phase. 

Dobrosavljevic, Miranda and coworkers [199-201] introduced the concept of an 
electronic Griffiths phase, a phase that occurs in strongly correlated electron systems 
in the vicinity of a disorder-driven metal insulator transition. These systems possess 
localized magnetic moments, that are either present from the outset (localized /- 
electrons in heavy fcrmion materials) or are self-generated due to the interplay between 
disorder and correlations. The electronic Griffiths phase is characterized by a power- 
law distribution of the Kondo temperatures Tr- of these moments. Integrating over 
this distribution leads to power-law singularities in many obscrvables including specific 
heat and magnetic susceptibility. More details can be found in the recent review [202] . 

Other rare region effects in quantum systems include the Lifshitz tails in the 
density of states in doped semiconductors [46,47]. They arc probably the oldest known 
rare region effects but not directly related to a phase transition. Analogous effects 
were also found in dirty superconductors [203-206]. Motrunich et al. [207] investigated 
Griffiths effects and quantum critical points in dirty one-dimensional superconductors; 
and Vafek et al. [208] discussed rare region effects at the flux-driven superconductor- 
metal transition in ultrasmall cylinders (see also reference [209]). Finally, Mildenberger 
et al. [210] identified a Griffiths phase in the thermal quantum Hall effect. 

6. Nonequilibrium phase transitions 

The nonequilibrium behavior of many-particle systems has attracted considerable 
attention in recent years. Of particular interest are continuous phase transitions 
between different nonequilibrium states. These transitions arc characterized by large 
scale fluctuations and collective behavior over large distances and times very similar to 
the behavior at equilibrium critical points. Examples of nonequilibrium transitions can 
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be found in population dynamics and epidemics, chemical reactions, growing surfaces, 
and in granular flow and traffic jams (for recent reviews see, e.g., references [21f 2f 7]). 
In this section we demonstrate that quenched spatial disorder at such nonequilibrium 
phase transitions leads to rare region effects that are analogous to those at equilibrium 
classical and quantum phase transitions. 

6.1. Contact process 

A prominent class of nonequilibrium phase transitions separates active fluctuating 
states from inac;tive. absorbing states where fluctuations cease entirely. Recently, much 
eflFort has been devoted to classifying possible universality classes of these absorbing 
state phase transitions [215,216]. The generic universality class is directed percolation 
(DP) [218]. According to a conjecture by Janssen and Grassberger [219,220], all 
absorbing state transitions with a scalar order parameter, short-range interactions, 
and no extra symmetries or conservation laws belong to this class. Examples include 
transitions in catalytic reactions [221], interface growth [222], or turbulence [223]. In 
the presence of conservation laws or additional symmetries, other universality classes 
can occur, e.g., the parity conserving class [224-226] or the Z2-symmetric directed 
percolation (DP2) class [227-230]. 

The contact process [231] is a prototypical system in the directed percolation 
universality class. It can be interpreted, e.g., as a model for the spreading of a disease. 
The contact process is defined on a d-dimensional hypercubic lattice. Each lattice site 
X can be active (occupied by a particle) or inactive (empty) . In the course of the time 
evolution, active sites can infect their neighbors, or they can spontaneously become 
inactive. Specifically, the dynamics is given by a continuous-time Markov process 
during which particles are created at empty sites at a rate An/ (2d) where n is the 
number of active nearest neighbor sites. Particles are annihilated at rate fi. The ratio 
of the rates A and fj, controls the behavior of the system. Without loss of generality 
we will set the annihilation rate /i to unity, unless noted otherwise. 

For small birth (infection) rate A, annihilation dominates, and the absorbing state 
without any particles is the only steady state (inactive phase). For large A, there is a 
steady state with finite particle density (active phase). The two phases are separated 
by a nonequilibrium phase transition in the directed percolation universality class at 
some critical value Aj?. The central quantity in the contact process is the average 
density of active sites at time t 



where nx(t) is the particle number at site x and time t, L is the linear system size, and 
(. . .) denotes the average over all realizations of the Markov process. The longtime 
limit of this density (i.e., the steady state density) 



is the order parameter of the nonequilibrium phase transition. 

Close to the nonequilibrium phase transition, the density shows scaling behavior 
very similar to that of equilibrium transitions discussed in section 2.2. Specifically, 
the order parameter pstat varies according to the power law 




(83) 
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Pstat ~ (A - Ac)'^ ~ r' 



(85) 
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where r = (A — \c)/\ is the dimensionless distance from the critical point, and /3 is 
the critical exponent of the particle density. In addition to the average density, we 
also need to characterize the length and time scales of the density fluctuations. Close 
to the transition, correlation length and correlation time diverge as 

^± ~ Irl-^- , ~ Ci, (86) 

i.e., the dynamical scaling is of power-law form with dynamical exponent 2;.|| 
Consequently the scaling form of the density as a function of r, the time t and the 
linear system size L reads 

p(r, t, L) = 6-'^/"^p(r6i/"^ , tb-\ Lb-') . (87) 

Here, b is an arbitrary dimensionless scaling factor. 

Two important quantities arise from initial conditions consisting of a single active 
site in an otherwise empty lattice. The survival probability Pg describes the probability 
that an active cluster survives when starting from such a single-site seed. For directed 
percolation, the survival probability scales exactly like the density 

Ps{r, t, L) = 6-^/"^P«(r6^/"^ , tb'^ Lb'') . (88) 

Thus, for directed percolation, the three critical exponents /?, z/^ and z completely 
characterize the critical point. ^ The pair connectedness function C(x',f',x, t) = 
(nx'(t') ^x(i)) describes the probability that site x' is active at time t' when starting 
from an initial condition with a single active site at x and time t. For a clean system, 
the pair connectedness is translationally invariant in space and time. Thus, it only 
depends on two arguments C{x,t' ,x,t) = C{x' — x,t — t'). Because C involves a 
product of two densities, its scale dimension is 2(3/v_\_, and the full scaling form reads 
(as long as hyperscaling is valid) 

C(r,x,i,L) = r2'3/''^C(r6l/'^-L,xfe-^^r^L5^l) . (89) 

The total number of particles N when starting from a single seed site can be obtained 
by integrating the pair connectedness function C over all space. This leads to the 
scaling form 

N{r,t,L) = b-^f^/^^+'^Nirb'/"^ ,tb-'' ,Lb-') . (90) 

At the critical point, r = 0, and in the thermodynamic limit, L — > oo, the 
above scaling relations lead to the following predictions for the time dependencies 
of observables: The density and the survival probability asymptotically decay like 
p{t) ~ t~^ and Ps{t) ~ t~^ with 5 = P/{i'±z). In contrast, the number of 
particles in a cluster starting from a single seed site increases like N{t) ~ t® where 
G = d/ z — 20/{p^z) is the so-called critical initial slip exponent. Highly precise 
estimates of the critical exponents for clean one-dimensional directed percolation have 
been obtained by series expansions [232]: (3 = 0.276486, vi_ = 1.096854, z = 1.580745, 
5 = 0.159464, and 6 = 0.313686. 

II Following the convention in the literature on the contact process, we denote the spatial correlation 
length and its critical exponent with a subscript ±. 

^ At more general absorbing state transitions, e.g., with several absorbing states, the survival 
probability scales with an exponent /3' which may be different from /3 (see, e.g., [215]) 
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6.2. Contact process with uncorrelated disorder 

Quenched spatial disorder can be introduced into the contact process by making the 
birth rate A or/and the death rate /i random functions of the lattice site x. In this 
subsection, we discuss point defects, i.e., spatially uncorrelated disorder. 

The investigation of disorder effects on spreading transitions in the directed 
percolation universality class actually has a long history, but a coherent picture has 
been slow to emerge. According to the Harris criterion [25], directed percolation is 
unstable against weak disorder because the clean spatial correlation length critical 
exponent z/^ violates the inequality diy± > 2 in all dimensions > 4. (The exponent 
values are iy± f« 1.097 (ID), 0.73 (2D), and 0.58 (3D) [215]). A field-theoretic 
renormalization group study [233] confirmed the instability of the directed percolation 
critical fixed point. Moreover, no new critical fixed point was found. Instead 
the renormalization group displays runaway flow towards large disorder, indicating 
unconventional behavior. Early Monte-Carlo simulations [234] showed significant 
changes in the critical exponents while later studies [235,236] of the two-dimensional 
contact process with dilution found logarithmically slow dynamics in violation of 
power-law scaling. In addition, slow dynamics was found in whole parameter regions 
in the vicinity of the phase transition [234, 237-240] . 

6.2.1. Griffiths singularities: Optimal fl,uctuation arguments. A first qualitative 
understanding of disorder and rare region effects in the disordered contact process 
can be gained by optimal fluctuation arguments similar to that in section 4.2. For 
definiteness, we assume the infection rate depends on the source site only and its 
probability distribution is binary, 



The impurity density p and their strength c arc constants between and 1. The 
impurities locally reduce the birth rate; the nonequilibrium transition will therefore 
occur at a value Ac that is larger than the clean critical birth rate A^. 

The inactive phase, A < Ac, can be divided into two regions. For birth rates 
below the clean critical point, A < A^, the behavior is conventional because non of 
the rare regions are locally in the active phase. The system approaches the absorbing 
state exponentially fast in time. The decay time increases with A and diverges as 
I A — Acl"^''^ where z and i'± are the exponents of the clean critical point [56,66]. In 
contrast, in the Griffiths region, i.e., for birth rates between the clean and the dirty 
critical points, Aj! < A < Ac, rare spatial regions devoid of impurities are locally in the 
active phase. Because they are of finite size, they cannot support a non-zero steady 
state density but their decay is very slow because it requires a rare, exceptionally large 
density fluctuation. 

Following Noest [234, 237] , the contribution of these rare regions to the time 
evolution of the total density can be estimated as follows. The probability w for finding 
a rare region of linear size -Lrr devoid of impurities is given by w;(Lrr) ^ cxp(— pL^p) 
where p = — ln(l — p) is a nonuniversal constant. The long-time decay of the density is 
dominated by these rare regions. To exponential accuracy, the rare region contribution 
to the density can be written as 



P[A(x)] = (1 - p) ^[A(x) - A] +p^[A(x) - cA] . 



(91) 




(92) 
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where t{L^b.) is the decay time of a rare region of size ^rr. Let us first discuss the 
behavior at the clean critical point, A^, i.e., at the boundary between the conventional 
inactive phase and the Griffiths region. At this point, the decay time of a single, 
impurity-free rare region of size Lrr scales as r(LRR) ~ L^^i ^ follows from finite 
size scaling [22]. Here, z is the clean dynamical exponent. A saddle-point analysis 
leads to a stretched exponential time dependence, 

Inp(t) - -p^-/(d+^-) t'^/(d+^) , (93) 

in contrast to the simple exponential decay as for A < A^. Inside the GrifHths region, 
i.e., for A° < A < Ac, the decay time of a single rare region depends exponentially on 
its volume t(Lrr) ^ exp(aLj|j^) because a coordinated fluctuation of the entire rare 
region is required to take it to the absorbing state [234,241]. The constant a vanishes 
at the clean critical point A° and increases with increasing A. Repeating the saddle 
point analysis of the integral (92) for this case, we obtain a power-law decay of the 
density 

pit) ~ t-p/" = r-^/^' (94) 

where z' = da/p is a customarily used nonuniversal dynamical exponent in the Griffiths 
region. Thus, the optimal fluctuation arguments predict strong power-law Griffiths 
effects for the disordered contact process. The reason is the exponential increase of 
the decay time r(LRR) with the size of the rare region that compensates for their 
exponential rarity. 



6.2.2. Strong- disorder renormalization group. An important step towards under- 
standing spatial disorder effects on the directed percolation transition has recently 
been made by Hooyberghs et al. [41,42]. These authors used the Hamiltonian formal- 
ism [242] to map the one-dimensional disordered contact process onto a random quan- 
tum spin chain. They then applied a version of the Ma-Dasgupta-Hu strong-disorder 
renormalization group [36] discussed in section 5.3 and showed that the transition is 
controlled by an infinite-randomness critical point, at least for sufficiently strong dis- 
order. The basic ingredients of this analysis can be explained by directly discussing 
the birth and death rates without reference to the equivalent quantum system. 

Let us consider a one-dimensional contact process where both the annihilation 
rate fii at site i and the infection rate Xij between neighboring sites i and j are 
random functions of the lattice sites with broad probability distributions. Following 
the strong-disorder renormalization group approach, we identify the largest rate in 
the system J7 = max(/j,j , A^ j ) . If the largest rate is an annihilation rate, say 1^2, site 
2 is almost always empty and can be decimated out. However, "virtual occupations" 
of site 2 generate a renormalized infection rate between sites 1 and 3 that can be 
estimated as follows: To propagate the infection from site 1 to site 3, first it has to 
propagate from site 1 to site 2 which happens with rate A12. However, site 3 can 
only be reached, if the infection propagates from site 2 to site 3 before site 2 decays. 
Since the probability distributions are broad, ^2 ^ A23, and the probability for this 
to happen is X23/H2. Thus, the renormalized infection rate between sites 1 and 3 is 

Al3 = Ai2A23/Ai2 • (95) 

If the largest rate is an infection rate, say A23, which connects sites 2 and 3, the two 
sites will (almost) always be in the same state, either both active or both inactive. 
They can therefore be replaced by a new single site with twice the weight (moment). 
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To estimate the renormalized annihilation rate, we note that total annihilation of 

the cluster can occur via two processes: (i) site 2 decays (with rate 112), then site 3 
decays before site 2 gets reinfected. Since A23 ^ ^3, the probability for this is /X3/A23. 
Thus the total rate for process (i) is /i2A*3/^23- Process (ii) consists of site 3 decaying 
first, followed by site 2 and makes an identical contribution to the rate. Thus, the 
renormalized death rate of the cluster is 

fl = 2/X2/^3/A23 . (96) 

Except for an extra factor of 2, the recursion relations (95) and (96) are 
completely equivalent to the recursion relations (62) and (63) in the strong 
disorder renormalization group analysis of the random transverse field Ising chain. 
Consequently, Hooyberghs et al. [41,42] also found equivalent results. The critical 
point in the disordered one-dimensional contact process is of infinite-randomness type, 
i.e., under renormalization, the relative with of the probability distributions of /Xj and 
Xij increase without limit. 

The scaling behavior at such an infinite randomness critical point is different from 
that of the clean contact process. Most importantly, the dynamics is extremely slow. 
The power-law scaling (86) gets replaced by activated dynamical scaling 

ln(6)~?t, (97) 

characterized by a new exponent tjj- This exponential relation between time and length 
scales implies that the dynamical exponent z is formally infinite. In contrast, the static 
scaling behavior remains of power law type. 

Moreover, at an infinite-randomness fixed point the probability distributions of 

observables become extremely broad, so that averages are dominated by rare events 
such as rare spatial regions with large infection rate. In such a situation, averages and 
typical values of a quantity do not necessarily agree. Nonetheless, the scaling form of 

the average density is obtained by simply replacing the power- law scaling combination 
tb~^ by the activated combination \n{t)b^'^' in the argument of the scaling function: 

p{r, ln(t), L) = b-^/'^^pirb^^"^ , ln(t)6-'^, Lfe-^) . (98) 

Analogously, the scaling forms of the average survival probability and the average 
number of sites in a cluster starting from a single site are 

Ps{r, Ht),L) = b-'^/'^^Psirb^'"^ , ln(i)6-'^, Lb'^) (99) 
N{rM{t),L) =b-'^P'''^+'^N{rb^'''^M{t)b-'l',Lb-^) . (100) 

These activated scaling forms lead to logarithmic time dependencies at the critical 
point (in the thermodynamic limit). The average density and the survival probability 
asymptotically decay like p{t) ^ [ln(t)]"^ and P,(f) [ln(t)]"^ with 5 = l3/{v±tp) 
while the average number of particles in a cluster starting from a single seed site 
increases like N{t) - [ln(i)]^ with 6 = d/ip - 2f3/{iyj_ip). Within the strong 
disorder renormalization group approach the critical exponents of the disordered one- 
dimensional contact process can be calculated exactly. The numerical values are 
/? = 0.38197, = 2, V = 0.5, S = 0.38197, and 6 = 1.2360. In spatial dimensions 
larger than one, the renormalization group can only be implemented numerically. As 
was discussed in section 5.4, this has been done for the strong-disorder renormalization 
group of the two-dimensional random transverse field Ising model [26]. Using the 
analogy between the renormalization group of the Ising model and the contact process, 
Hooyberghs et al. [42] therefore concluded that the disordered contact process in two 
dimensions will also show an infinite randomness critical point. 
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Figure 11. Left:p~^/'' vs. ln{t) for a system of 10^ sites witli p = 0.3 and c = 0.2. 
The data are averages over 480 runs, eacli with a different disorder realization. 
The filled circles mark the critical birth rate Ac = 5.24, and the straight line is a 
fit of the long-time behavior to eq. (101). Right: Log-log plot of the density time 
evolution in the Griffiths region for systems with p = 0.3, c = 0.2 and several birth 
rates A. The system sizes are lO'^ sites for A = 3.5, 3.7 and 10® sites for the other 
A values. The straight lines are fits to the power law p(t) ~ t~^/^ predicted in 
eq. (94). Inset: Dynamical exponent z' vs. birth rate A (from [43]). 



In addition to the infinite-randomness behavior at the critical point, the strong- 
disorder renormalization group also yields strong power-law rare region effects in the 
Griffiths phase in agreement with the optimal fluctuation arguments presented in 
section 6.2.1. 

6.2.3. Computer simulations. Early Monte-Carlo simulations of the disordered 
contact process [234-237] gave somewhat inconclusive results ranging from a 
modification of the critical exponents to logarithmical violations of scaling. This 
was partially due to the fact that a consistent theory was not available at that 
time. Hooyberghs et al. [42] performed both Monte-Carlo and density matrix 
renormalization group (DMRG) studies of rather small systems {L < 24). They 
concluded that the phase transition is indeed controlled by the infinite randomness 
fixed point discussed in section 6.2.2, but only if the initial (bare) disorder is large 
enough. For weaker initial disorder they predict noniiniversal continuously varying 
exponents with either power-law or activated dynamical scaling. 

Recently, Vojta and Dickison [43] have performed a detailed large-scale Monte- 
Carlo study of the one-dimensional contact process with binary disorder distribution 
(91). They studied system sizes up to L = 10*" and very long-times up to f = 10^ 
which is three orders of magnitude in t longer than previous simulations. A first 
set of calculations started from a full lattice and followed the time evolution of the 
average density. To test the strong-disorder renormalization group theory, the time 
dependence of the density can be plotted according to the predicted activated scaling 
law 

Pit) ^ Mt)]-', (101) 

with 6 = 0.38197. The left panel of figure 11 shows the resulting graph for a system of 
10^ sites with p = 0.3 and c = 0.2. The evolution of the density at the critical birthrate 
A = 5.24 follows eq. (101) over almost six orders of magnitude in time. Therefore, the 
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Pc P' 

Figure 12. Schematic phase diagram of a site diluted contact process as function 
of impurity concentration p and birth rate A. There is a muhicritical point at 
p = Pc and A = At (from [243]). 

critical point of the disordered contact process is indeed of infinite-randomness type. 
Analogous calculation were performed for two other sets of parameters (p = 0.3 with 
c varying from 0.2 to 0.8 as well as c = 0.2 with p varying from 0.2 to 0.5). For all 
parameters (including the cases of weak disorder), the critical point is characterized by 
the logarithmic density decay (101) with a universal exponent 5 = 0.38197. In addition 
to the critical point, Vojta and Dickison [43] also studied the Griffiths region between 
the clean critical birthrate, A° — 3.298 and the dirty critical birthrate Ac = 5.24. The 
right panel of figure 11 shows a double-logarithmic plot of the density time evolution 
for birth rates A — 3.5 . . . 5.1 and p — 0.3, c = 0.2. For all birth rates A shown, p{t) 
follows (94) over several orders of magnitude in p (except for the largest A where 
we could observe the power law only over a smaller range in p because the decay is 
too slow). The nonuniversal dynamical exponent z' can be obtained by fitting the 
long-time asymptotics of the curves to eq. (94). The inset of the figure shows z' as a 
function of the birth rate A. As predicted, z' increases with increasing A throughout 
the Grifhths region with an apparent divergence around A = Ac = 5.24. 

These large-scale simulations of the one-dimensional contact process with point 
defects thus provide strong evidence that the critical point is of infinite-randomness 
type with universal critical exponents (even for weak bare disorder) . The critical point 
is accompanied by strong Griffiths singularities characterized by non-universal power- 
law decay of the density. These results are in excellent agreement with the general 
classification of dirty phase transitions suggested in section 3.3. 

6.2.4- Contact process on a randomly diluted lattice. Quenched spatial disorder can 
also be introduced into the contact process by site or bond dilution of the underlying 
lattice. In space dimensions d > 2, this leads to the phase diagram shown in figure 
12 which is very similar to that of a diluted quantum magnet discussed in section 5.4. 
(In contrast, in d = 1, any dilution immediately destroys the active phase.) For small 
impurity concentrations below the percolation threshold of the lattice, p < Pc, the 
active phase survives, but the critical birth rate increases with p (to compensate for 
the missing neighbors). Right at the percolation threshold the active phase survives 
on the infinite percolation cluster for A > A*. For p > pc, no active phase can exist 
because the lattice consists of disconnected clusters of finite size that do not support 
a steady state density of active sites. 
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The contact process on a site-diluted lattice therefore has two nonequilibrium 
phase transitions, separated by a multicritical point at (pc,A»). For p < Pc, the 
transition (marked by "a" in Fig. 12) is expected to be in the universality class of 
the generic disordered contact process discussed in the last subsections. In contrast, 
the phase transition "b" across the percolation threshold of the lattice is expected to 
show different critical behavior. 

Vojta and Lee [243] showed that the interplay of geometric criticality and dynamic 
fluctuations of the contact process does indeed lead to a novel universality class for 
this nonequilibrium phase transition. Even though the transition is driven entirely by 
the geometry of the lattice, the dynamical fluctuations of the contact process enhance 
the singularities in all quantities involving dynamic correlations. By combining results 
from classical percolation theory [147] with the finite-size properties of a supercritical 
contact process, Vojta and Lee [243] determined the critical behavior analytically. 
They found activated dynamical scaling In ~ with the critical exponent ip being 
equal to the fractal dimension of the critical percolation cluster, ip = Df. As a result, 
the critical behavior is very similar to that of the infinite-randomness critical point 
discussed above. Specifically, the long-time decay of the density p of active sites at 
p = Pc is logarithmically slow, p(t) ~ [\n{t/to)]~^ . The exponent 5 = (idi^cDf) is 
determined by D j together with the order parameter and correlation length exponents, 
f3c and Vc, of the lattice percolation transition. In contrast to the enhanced dynamical 
singularities, the exponents of static quantities like the steady state density Pst and 
the spatial correlation length are identical the corresponding lattice percolation 
exponents, Psiip) ~ I'i"' and ^ [rj'^" where r =p — Pc measures the distance from 
the percolation threshold. Off criticality, i.e., away from the percolation threshold, 
there are strong rare region effects characterized by a non-exponential density decay 

Pit) ~ (t/to)-'^/^' for p> Pc and p{t) - pst ^ g-I^'*/^"' for p < pc. The 
nonuniversal exponents z' and z" diverge as for p ^ Pc- 

6.3. Contact process with extended defects 

So far, our discussion of disorder effects on the phase transition in the contact process 
has been restricted to uncorrelated or short-range correlated disorder (i.e., point 
defects). In this subsection we consider spatially extended (linear or planar) defects. 
Prom the general arguments in section 3.2, we expect the disorder correlations to 
enhance the impurity effects. 

Vojta [65] studied a contact process with quenched spatial disorder perfectly 
correlated in c?cor > space dimensions but uncorrelated in the remaining c?ran = 
d — dcor space dimensions. In such a system, the rare regions arc infinite in licor 
dimensions but finite in dran dimensions. This is a crucial difference from systems 
with uncorrelated disorder, where the rare regions are finite. Because the contact 
process has a phase transition in all dimensions larger than zero, an infinite rare region 
can undergo a real phase transition independently of the rest of the system. Those 
rare regions that are locally in the ordered phase will have a true nonzero stationary 
density, even if the bulk system is still in the inactive phase. The resiilting global 
phase transition is thus smeared by the same mechanism that smears the equilibrium 
transition in an Ising model with planar defects discussed in section 4.6 or the itinerant 
Ising magnet of section 5.6. As the birth rate A is increased, the order parameter 
develops very inhomogeneously in space with different parts of the system becoming 
active independently at different values of A. 
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The leading behavior in the tail of the smeared transition can again be determined 
by optimal fluctuation arguments [65]. Here wc summarize the results for a binary 
probability distribution P[A(Xran)] = (1 — p) <5[A(Xran) — M +P^[H^Ta.n) — cA] with 
constants p and c between zero and one. Xran is the projection of the position vector x 
on the uncorrelated directions. The probability w for finding a rare region of linear size 
-C'RR devoid of impurities is, up to pre-exponential factors, given by ui ~ exp(— pL^Jg") 
with p = — ln(l — p). According to finite-size scaling [22], such a region undergoes 
a true phase transition to the active phase at Xc{Lkr) = AJ? + AL^^ where <j> is the 
clean (d-dimensional) finite-size scaling shift exponent and ^ is a constant. If the 
total dimensionality d = dcor + c^ran < 4, hyperscaling is valid and (j) = l/i'± which we 
assume from now on. The total (average) density p at a certain A > A° is obtained by 
summing over all active rare regions. In this way, the stationary density develops an 
exponential tail, 

p(A)~exp(-B(A-AO)-*-''-) , (102) 

reaching to the clean critical point A^. Here, B = pA'^'^""-^ is a constant. Analogous 
arguments can be made for the survival probability P{\) of a single seed site. If the 
seed site is on an active rare region it will survive with a probability that depends on 
A with a power law. If is not on an active rare region, the seed will die. To exponential 
accuracy the survival probability is thus also given by (102). 

We now turn to the dynamics in the tail of the smeared transition. The long- 
time behavior of the density is dominated by the slow decay of the rare regions while 
the bulk contribution decays exponentially. Proceeding in analogy to section 4.8, 
Vojta [65] showed that the leading long-time decay of the density at the clean critical 
point, A = A^, follows a stretched exponential, 

where z is the clean dynamical exponent of the rf-dimensional bulk system. For A > A°, 
the behavior changes. At times below a crossover time ta; ~ (A — A°)~('^'a"+^)''-L the 
density decay is still given by the stretched exponential (103). For times larger than tx 
the system realizes that some of the rare regions are in the active phase and contribute 
to a finite steady state density. The approach of the average density to this steady 
state value is characterized by a nonuniversal power-law. 

p{t) - p((X)) - t-^ . (104) 

To test the predictions of the optimal fluctuation theory, Dickison and Vojta [66] 
have performed large-scale Monte-Carlo simulations of a contact process on a square 
lattice. The disorder consists of linear defects, i.e., the local birthrate A(x) depends 
only on one of the coordinates and is perfectly correlated in the other. The disorder 
distribution is binary. 

The left panel of figure 13 shows a comparison of the stationary density Pst as 
a function of A between the clean system and a dirty system with p = 0.2, c = 0.2. 
The clean system (p = 0) has a sharp phase transition with a power-law singularity 
of the density, pst ~ (A — A°)^ at the clean critical point A° « 1.65 with /3 0.58 
in agreement with the literature [235]. In contrast, in the dirty system, the density 
increases much more slowly with A after crossing the clean critical point. This suggests 
either a critical point with a very large exponent /? or exponential behavior. To test 
the prediction (102) of the optimal fluctuation theory, the right panel of figure 13 
shows In as a function of (A — A^)""^-^ for several impurity concentrations p. (The 
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Figure 13. Left: Stationary density pst as a function of birtli rate X for a 
clean system and a system witli impurity concentration p = 0.2 and strength 
c = 0.2. The data are averages over 200 disorder realizations and the system 
size is L = 1000. Right: Logarithm of the stationary density /Ogt as a function of 
(A — A")^'-'^ = (A — A")^" ''^* for several impurity concentrations p and L = 3000. 
The straight lines are fits to eq. (102) (from [66]). 
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Figure 14. Logarithm of the density at the clean critical point Aj! as a function of 
^i/(i+z) _ ^0.362 £qj, ^ _ 2 and several impurity concentrations {p = 0.2, . . . ,0.4 
from top to bottom) and system size L = 3000. The long-time behavior follows a 
stretched exponential Inp = —Et^-^^'^ (from [66]). 



clean two-dimensional spatial correlation length exponent is v±_ = 0.734 [244].) The 
data show that the density tail is indeed a stretched exponential, following the law 
(102) over at least two orders of magnitude in Fits of the data to eq. (102) can be 
used to determine the decay constants B. As predicted by extremal statistics theory, 
the decay constants depend linearly on p — — ln(l — p). 

In addition to the stationary density, Dickison and Vojta [66] also studied its time 
evolution in the tail of the smeared transition. Figure 14 shows the time dependence 
of the density right at the clean critical point. The data follow the predicted stretched 
exponential behavior (103) with z = 1.76 [244] over more than three orders of 
magnitude in p. Additional calculations have been performed for A > A°. In this 
parameter region, the stationary density is nonzero. The approach of the density to 
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this nonzero stationary value follows a nonuniversal power-law as predicted in (104). 

6.4- Other nonequilibrium systems 

While the influence of quenched spatial disorder on the directed percolation 
universality class has been studied in great detail, other dirty absorbing state 
transitions have received much less attention. Hooyberghs et al. [42] studied disorder 
efi^ects on the transition of a generalized contact process with two absorbing states 
[230]. This model is in the DP2 [227] or parity conserving [224] universality class 
(the two coincide in one space dimension [215]). According to the Harris criterion, 
the clean transition is unstable with respect to weak disorder because the clean 
correlation length exponent is p± « 1.83 (see, e.g., [215]) violating the inequality 
di'± > 2. Hooyberghs ct al. applied the strong-disorder renormalization group to a 
weakly disordered generalized contact process with quenched spatial disorder. They 
found that the renormalization group scheme does not work; the flow is not towards 
larger disorder. They concluded that the transition is unlikely to be controlled by 
an infinite disorder fixed point. Very recently. Odor and Menyhard [245] reported 
numerical simulations of a cellular automaton in the parity conserving universality 
class. They found that weak disorder does not change the critical exponents from 
their clean values. How this can be reconciled with the violation of the Harris criterion 
remains a question for future research. 

Let us also briefly mention another important nonequilibrium model displaying 
rare region efl^ects, even though it is not directly related to an order-disorder phase 
transition, viz. Sinai's model of diffusion with random local bias [246]. In the absence 
of a bias, the mean-square displacement in this model grows logarithmically slowly 
in time, while for small bias it grows with a nonuniversal power-law. Igloi, Turban 
and Rieger [119,247] pointed out a connection between the Ising model and directed 
walks. Later, Igloi and Rieger [248] established an exact relation between the statistical 
properties of anomalous diffusion and the properties of random quantum spin chains. 
At about the same time Fisher, Le Doussal and Monthus used a version of the strong- 
disorder renormalization group to determine the exact asymptotic long-time behavior 
of the Sinay model [138, 139]. These papers started a stream of work on this topic. 
A detailed discussion is beyond the scape of this review, for more details see, e.g., 
reference [121]. 

7. Discussion and conclusions 

In this section we summarize the main theoretical ideas presented in the review. To 
this end, we organize the examples presented in sections 4, 5, and 6 according to 
the proposed classiflcation of rare region effects, and we relate our results to those 
of other methods. We also discuss the generality of the rare region classification, its 
limitations, the experimental situation, and some related questions. 

7.1. Summary: Classification of rare region effects 

This topical review has been devoted to rare region effects at order-disorder transitions 
between conventional phases. In section 3.3 we have introduced a classification of these 
effects according to the effective dimensionality c^rr of the rare regions [67] . Its basic 
idea is as follows: The probability w for finding a rare region of volume Vr,r (in the 
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random directions) decreases exponentially with its size, w ~ exp(— ^Vrr). This must 

be compared to how rapidly the contribution of a single region to thermodynamic 
quantities increases with its size. As was explained in section 3.3, three cases can be 
distinguished. 

Class A: If the effective dimensionality of the rare regions is below the lower 
critical dimension of the problem, dnn < d~ , the contribution of a rare region to 
thermodynamic quantities only increases like a power of its size. Therefore, Griffiths 
effects arc expected to be exponentially weak and the dirty critical point will display 
power-law scaling. 

Class B: For rfRR = d~ , rare region contributions increase exponentially with the 
size of the rare region resulting in strong power- law Griffiths singularities. The dirty 
critical point is expected to be of infinite randomness type with activated scaling. 

Class C: If c^rr > d~ , an isolated rare region can undergo the phase transition 
independent of bulk and develop a static order parameter. This completely destroys 
the global phase transition by smearing. 

In the following we relate the explicit results discussed in sections 4, 5, and 6 
to this classification. We start with the classical systems of section 4. The classical 
randomly diluted Ising model introduced in section 4.1 as well as the classical random- 
Tc Landau-Ginzburg- Wilson theory of section 4.3 have spatially uncorrelated disorder. 
The rare regions arc of finite size in all directions and their effective dimensionality 
is cJrr = 0. Since the lower critical dimension of the Ising model is = 1, 
these systems fall into class A. In agreement with the classification, their Griffiths 
singularities arc exponentially weak, as was discussed in sections 4.2 and 4.3. In 
contrast, the properties of the dirty critical point are still controversial. On the 
one hand, high precision numerical simulations (see, e.g., [31,249]) yield conventional 
power-law critical behavior and no indications of exotic properties, in agreement with 
early perturbative renormalization results [250,251]. Moreover, the exponent values 
are in quantitative agreement with re-summed high-order perturbative renormalization 
group calculations [252- 256]. On the other hand, Dotsenko et al. [257,258] analyzed 
the stability of the perturbative fixed points with respect to replica symmetry breaking 
and found them unstable, suggesting unconventional behavior. However, two-loop field 
theoretical studies [259, 260] did not find such instabilities with respect to replica- 
symmetry breaking. Clearly, more work is necessary to determine whether the critical 
point in the Ising model with uncorrelated disorder indeed follows all of the predictions 
of class A. 

If the disorder is perfectly correlated in (icor dimensions, the rare regions are 
infinite in the corresponding directions but finite in the remaining dran = d — dcor 
directions. Thus, there effective dimensionality is dpR = dcor- The McCoy-Wu 
model discussed in section 4.5, a two-dimensional Ising model with linear defects, 
falls into class B because c^rr = d~ = 1. In agreement with the classification, this 
model displays strong power-law Griffiths singularities and an infinite-randomness 
critical point. Using Monte-Carlo simulations [40] and a numerical implementation of 
the strong-disorder renormalization group [26], the same behavior was also found in 
two dimensions. As was discussed in section 4.5, perturbative renormalization group 
calculations [106,107,261] cannot capture the nonperturbative rare region physics, they 
predict conventional behavior in disagreement with the exact results. In this context, 
an investigation of replica-symmetry breaking in systems with extended defects would 
be of great interest. 

In an Ising model with planar defects, the rare regions have dimension c^rr = 2. 
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This is above the lower critical dimension, and the model falls into class C. The 
cxplicdt results discussed in sciction 4.6 show that the global phase transition is indeed 
smeared [63] as was also confirmed by large-scale Monte Carlo simulations [64]. 

In the case of zero-temperature quantum phase transitions, the imaginary time 
direction has to be included when calculating the effective rare region dimension rfrm- 
Therefore, the rare regions in the random transverse field Ising model introduced 
in section 5.1 are one-dimensional, dRR = 1, and the model falls into class B. The 
strong disorder rcnormalization group results [26,34,35] and Monte-Carlo simulations 
[39,40,116] discussed in sections 5.1 to 5.4 are in complete agreement with this 
classification: The critical point is of infinite-randomness type with activated scaling. 
It is accompanied by strong power-law Griffiths singularities. 

At the quantum phase transition of the diluted bilayer Heisenberg quantum 
antiferromagnet discussed in section 5.5, the effective dimensionality of the rare regions 
in rfRp = 1, too. However, the lower critical dimension of the Heisenberg universality 
class isd~ =2. Therefore, this system falls into class A. This agrees with the results of 
large-scale Monte-Carlo simulations [70] that show a conventional critical point with 
power-law scaling and universal exponents. 

In systems with dissipative order parameter dynamics like the itinerant magnets 
of section 5.6, the behavior gets modified because the effective interaction in imaginary 
time direction takes a long-ranged l/r^ form. The rare regions therefore correspond to 
one- dimensional classical systems systems with 1/x^ interaction. For Ising symmetry, 
the rare regions can undergo the phase transition independently from the bulk, and 
the system is in class C. This results in a smeared global phase transition [60]. One- 
dimensional Heisenberg models with 1/x^ interaction are exactly at their lower critical 
dimension. Consequently, disordered itinerant Heisenberg magnets fall into class B 
and show power-law quantum Griffiths effects [67] (at least as long as long-range 
interactions between the rare regions can be neglected [73]). As is to be expected, a 
straight-forward perturbative analysis [179] of the transition in a disordered itinerant 
magnet does not predict any unconventional properties. However, Narayanan et al. 
[262,263] included the rare regions in a perturbation renormalization group calculation 
in an approximate way. They found that the rare regions destabilize the conventional 
fixed points but they could not resolve the ultimate fate of the transition. 

In section 6, we have demonstrated that the rare region classification also applies 
to nonequilibrium phase transitions in the directed percolation universality class. In 
the contact process with uncorrclatcd disorder (point defects) , the effective dimension 
of the rare regions is cIrr = 1 including the time dimension. This is identical to the 
lower critical dimension of directed percolation, and thus the system is in class B. In 
agreement with the classification, the strong-disorder rcnormalization group as well as 
Monte-Carlo simulations yield an infinite-randomness critical point and strong power- 
law Griffiths effects. If the defects are extended, the rare region dimension is c^rr > 1 
placing the system in class C. The classification thus predicts a smeared transition in 
agreement with the explicit results in section 6.3. 

In summary, all results presented in sections 4, 5, and 6 fit into the general 
classification of rare region effects at weakly disordered order-disorder transitions, as 
introduced in section 3.3. 
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7.2. Experiments 

This subsection is devoted to experimental observations of the rare region effects 
discussed in this review. 

7.2.1. Classical phase transitions. As explained in section 4, the thermodynamic 
rare region effects at generic classical (i.e., thermal) phase transitions are exponentially 
weak. As was first argued by Imry [52] they are probably unobservable in experiments, 
because the probability for finding sufficiently large rare regions is extremely small even 
in a macroscopic sample of about 10^^ particles. Indeed, to the best of our knowledge, 
the thermodynamic Griffiths singularities derived in sections 4.1 to 4.3 have never 
been experimentally verified. 

In contrast to the thermodynamics, the long-time dynamical behavior of generic 
classical systems is predicted to be dominated by the rare regions, see section 4.4. 
Lloyd and Mitchell [264] investigated the dynamics above the Neel temperature 

= 210K of the dilute Heisenberg antiferromagnet KMno.3Nio.7F3 using inelastic 
neutron scattering. Below the transition temperature = 246K of pure KNiFs they 
found very slow relaxation compatible with a long-time tail of stretched exponential 
form. However, the data were not sufficient for a quantitative confirmation of (36). 

In recent years, there have been a number of systems in which experiments have 
revealed unusual cluster or domain structures. Many of these phenomena have been 
related related to Griffiths phases, at least at a qualitative level. Binek, Kleemann, and 
coworkers [265-267] found domain-like antiferromagnetic correlations in the Ising-type 
antiferromagnets FeCl2 and FeBr2 when exposed to an axial field H at temperatures 
below the Neel temperature on zero field. They attributed theses structures to 
fluctuating distributions of demagnetizing fields and called the corresponding phase 
a field-induced Griffiths phase. However, it should be emphasized that the character 
of this domain phase is fundamentally different from the Griffiths phases discussed 
in this review. The latter are caused by static disorder fluctuations while the 
demagnetizing fields fiuctuates in time. Later, Binek and Kleemann also studied 
the diluted ferromagnet K2Cui_a;Zna;F4 and the diluted antiferromagnet Fei_a;Zna;F4 
[268] by measuring the magnetic susceptibility. They found deviations from the Curie- 
Weiss behavior of the corresponding clean systems, but no relation to the exponential 
Griffiths singularities discussed in sections 4.1 to 4.3 was established. Similar results 
were also reported for Rb2Coi_xMga;F4 [269]. CoUa ct al. [270] studied the dielectric 
response of a disordered antiferroelectric, Pbo.97Lao.o3(Zro.6oSno. 30 Tio. 10)03. They 
found nonanalytic contributions compatible with a Griffiths scenario in the nonlinear 
dynamic susceptibility but not in the dependence on dc fields. 

An important class of materials whose properties have been related to rare region 
effects are the manganites. They are members of a broader class of compounds where 
the electronic correlations play an important role. The manganites have rich phase 
diagrams exhibiting a variety of phases with unusual spin, charge, lattice and orbital 
order; for details see, e.g., the reviews [271,272] and the references therein. They 
display many exotic properties, most notably unexpectedly large magnetotransport 
properties. Moreover, even in the best available crystals, manganites are intrinsically 
inhomogeneous, a phenomenon sometimes called nanoscale phase separation. Salamon 
and coworkers [273, 274] proposed a connection between the unusual properties of 
doped LaMnOs and the Griffiths phase that arises when disorder suppresses the 
magnetic transition (see also reference [275] for a theoretical account). They found that 
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the magnetic susceptibility in Lao.TCao.sMnOs strongly increases before the system 
reaches the Curie temperature. They interpreted this upturn as a Griffiths singularity 
and analyzed it using the susceptibility matrix approach of Bray [55]. However, an 
explicit verification of the exponential classical Griffiths singularity was not found. 
Hartingcr et al. [276] and Deisenhofer et al. [277] found inhomogencous states in 
Lai_a;Sra;Mn03 that they attributed to Griffiths-like phases. The latter paper also 
pointed out that the disorder in the manganites is likely to be correlated. Because 
these disorder correlations enhance the rare region effects they may be responsible for 
making the Griffiths singularities in the manganites accessible to experiment. 

A Griffiths phase has also been predicted to occur in diluted magnetic 
semiconductors above the ferromagnetic transition temperature [278]. Recent 
measurements on MuxGei-x [279] are compatible with a scenario in which 
ferromagnetic long-range order occurs below some Tc and short-range order occurs 
between Tc and T*. (For x = 0.05, Tc w 12K and T* w 112K). However, the explicit 
form of the Griffiths singularities was not studied. 

7.2.2. Quantum phase transitions. While thermodynamic Griffiths singularities at 
generic classical transitions are very weak, looking for rare region effects at zero- 
temperature quantum phase transitions is a-priori more promising, because quantum 
Griffiths effects arc gencrically stronger than classical ones. 

When the prototypical transverse-field Ising magnet LiHoF4 is diluted by 
randomly replacing magnetic Ho^"*" with nonmagnetic Y^+, the long-range dipolar 
character of the interaction leads to frustration. In a transverse magnetic field, the 
result is therefore a quantum spin glass rather than a random quantum ferromagnet. 
LiHoi-2;Y2;F4 has been studied extensively since the seminal work of Wu et al. 
[114,115], but its complex properties are beyond the scope of this review. 

In recent years, there has been a lot of interest in the possibility of quantum 
Griffiths effects in the so-called heavy- fermion materials. These are intermetallic 
compounds of rare earth or actinide elements such as Cerium, Ytterbium or Uranium. 
The localized magnetic moment of the /-electrons of these atoms hybridizes with the 
delocalized conduction electrons generating quasi-particles with effective masses of 
100 to 1000 electron masses. Many of these systems can undergo magnetic quantum 
phase transitions as a function of pressure or doping or magnetic field. Moreover, many 
show non- Fermi liquid behavior, i.e., anomalous temperature dependencies of specific 
heat, magnetic susceptibility and resistivity, as was first shown for Y0.8U0.2Pd3 in 
1991 [280]. Since then, similar behavior has been found in a huge number of systems 
(for a detailed experimental review, see reference [281]). 

Castro Ncto, Castilla, and Jones [183] suggested that the anomalous temperature 
dependencies are a manifestation of power-law quantum Griffiths effects in the vicinity 
of a magnetic quantum phase transition. However, this early version of the theory 
neglected the damping of the magnetic modes due to the itinerant electrons. As was 
discussed in section 5.6 of this review, this damping qualitatively changes the low- 
temperature behavior. For Ising symmetry, which is relevant for the heavy-Fermion 
systcinis as most of them have sizeable spin anisotropics, the damping causes the 
rare regions to freeze leading to a smeared transition rather than quantum Griffiths 
effects [60-62]. Later, Castro Neto and Jones refined their theory by including the 
damping [282,283]. They found that for sufficiently weak damping there will be 
a crossover temperature T*. Power-law quantum Griffiths behavior can occur in a 
transient temperature window above T* but below a microscopic cutoff scale (like 
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the Fermi temperature). For the experimentally important example of the heavy 
fcrmion materials, different approaches disagree concerning the value of this crossover 
temperature. Castro Neto and Jones [282-284] suggest that the damping of the locally 
ordered rare regions is weak. Therefore, quantum Griffiths behavior should occur in 
a broad temperature range and can be the source of the non-Fermi liquid behavior 
observed in at least some of the heavy-Fermion systems. In contrast, Millis, Morr and 
Schmalian [62, 285] argue that the carrier-spin coupling in heavy fermions is large, 
restricting quantum Griffiths effects to a narrow temperature window (if any). A 
complete resolution of this question will likely come from experiments that give more 
direct access the the strength of the damping term. 

Recently, Johnston ct al. [286] observed slow non-exponential nuclear magnetic 
relaxation in NMR experiments on the d-electron heavy fermion system LiV204 with 
a small amount of magnetic defects. One possible explanation for this observation 
consists in the defects inducing magnetic order in large but finite volumes around the 
defect site. Johnston et al. [286] described these magnetic clusters by Griffiths physics. 

Let us also briefly mention some experiments on random quantum spin chains. 
Early examples include the organic crystals NMP-TCNQ [287], Qn(TCNQ)2 [288], 
and copper-doped TMMC [289]. The former two systems show low-temperature 
susceptibilities proportional T~" with a < 1 that have been interpreted as evidence for 
rare region physics. The latter system, TMMC, realizes an almost classical spin chain, 
presumably due to the large value of the spin involved {S = 5/2 for Mn^+). More 
recently, Payen et al. [290] reported susceptibility and magnetization for the doped 
S = 1 Haldanc spin system Y2BaNii_2;Zn2,05. Below 4K they observed power-law 
temperature dependencies with non-universal exponents suggesting a gapless phase 
due to the disorder. Power-law temperature dependencies have also been reported for 
the S = 5/2 spin chain system MnMgB205 above a spin glass freezing temperature of 
about 600mK [291]. Masuda et al. [292] reported inelastic neutron scattering on the 
S = 1/2 spin chain compound BaCu2(Sio.5Geo. 5)207. They found excellent agreement 
with the results of a strong-disorder renormalization group [207]. 

7.2.3. Nonequilibrium phase transitions. For nonequilibrium phase transitions, the 
experimental situation is even less satisfactory. The (clean) directed percolation 
universality class is ubiquitous in theory and simulations, but experimental 
verifications are strangely lacking [293]. To the best of our knowledge, the only 
(partial) verification so far has been found in the spatio-temporal intermittency in 
ferrofluidic spikes [294] . While disorder is often suggested as a reason for the striking 
absence of clean directed percolation scaling in at least some of the experiments, a 
systematic search for the infinite-randomness behavior discussed in section 6.2 and 
the accompanying power law Griffiths effects has not been performed. 

In summary, even though rare region effects have been invoked in the analysis of 
a variety of experiments, clearcut experimental verifications of Griffiths singularities 
discussed in sections 4, 5, and 6 have not yet been achieved. 

7.3. Conclusions 

Quenched spatial disorder can have dramatic consequences for the properties of 
continuous phase transitions and critical points. An important role is played by rare 
strong disorder fluctuations and the rare spatial regions where they occur. This topical 
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review has focused on transitions of the order-disorder type at which the quenched 
disorder docs not qualitatively change the properties of the phases but locally modifies 
the tendency towards one or the other. For this random- or random-mass type 
disorder (arguably the simplest type of quenched disorder), we have set up a general 
framework for understanding rare region effects in the vicinity of the "dirty phase 
transitions". According to the classification introduced in section 3.3, the overall 
phenomenology of the transition is determined by the effective dimensionality of the 
rare regions. 

Let us briefly point out the limitations of this approach. Even for order-disorder 
transitions between conventional phases, the effects of long-range spatial interactions 
arc not understood. First, the very concept of a rare region may require modifications 
when the order parameter of a locally ordered island generically develops power- 
law tails. Second, long-range interactions can significantly increase the coupling 
between the rare regions. As a result, the rare region effects in systems with long- 
range interactions are likely to be more complex than those discussed here. This is 
particularly important in metallic systems, were the presence of electronic soft modes 
often self-generates long-range order parameter interactions [14]. For instance, it has 
been suggested [73] that the long-range RKKY interaction in metals further enhances 
the rare region effects in itinerant Heisenberg magnets. 

Often, the presence of quenched disorder leads to the formation of new phases, 
such as the spin glass phases in classical magnets or the random-singlet states of certain 
antiferromagnetic spin chains, as briefly mentioned in section 5.7. These states, and 
phase transitions between them, are beyond the rare region classification in its present 
form. 

In recent years, there has been strong interest in "exotic" phases and phase 
transitions that do not follow Landau's order parameter paradigm (for a recent 
pedagogical overview, see, e.g., [295]. The influence of disorder on such transitions 
has not been studied in any detail, and it is likely very different from the behavior 
discussed here. 

Finally, Griffiths singularities are normally associated with continuous phase 
transitions, and we have also considered this case. Concerning first-order phase 
transitions, it has been known for a long-time that quenched disorder generally 
weakens the first-order nature. Building on earlier work by Imry and Wortis [296], 
Aizenman and Wehr [297] proved, that in dimensions d < 2, a. first-order transition 
turns into a continuous one upon introduction of arbitrarily weak (random-mass type) 
disorder. However, for d > 2, first order transitions can persist in the presence of 
quenched disorder. Not much is known about pretransitional phenomena at such 
first-order transitions (see, however, reference [298]). Exploring this in more detail 
would also be important from an experimental point of view, as some of the transition 
mentioned in the last section (e.g., in the manganites) appear to be of first order. 
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